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New challenge for bionics—brain-inspired computing 


Shan YU 


By definition, bionics is the application of biological mechanisms 
found in nature to artificial systems in order to achieve specific 
functional goals. Successful examples range from Velcro, the 
touch fastener inspired by the hooks of burrs, to self-cleaning 
material, inspired by the surface of the lotus leaf. Recently, a 
new trend in bionics—Brain-Inspired Computing (BIC)—has 
captured increasing attention. Instead of learning from burrs 
and leaves, BIC aims to understand the brain and then utilize its 
operating principles to achieve powerful and efficient 
information processing. 

In the past few decades, we have witnessed dramatic 
progress in information technology. Moore’s law, which states 
that transistor density in processors doubles every two years, 
has been proven true for the last 50 years. As a result, we now 
have miniature processors in small devices (e.g., phones) that, 
in terms of numerical calculation and memory storage, easily 
dwarf the brightest human mind. Given such a condition, which 
aspects of the brain can still enlighten us? 

First, we need more energy-efficient processors. Nowadays, 
supercomputers or large data centers contain thousands of 
cores/processors, with the energy consumption rate at the 
megawatt scale. This severely limits the use of computing 
power in embedded (e.g., small, smart devices) and long 
distance (e.g., Mars rover) applications. In addition, with further 
extrapolation of Moore's law, the energy density of a 
microprocessor will become so high that it will start to melt. In 
fact, this is an important reason why it is believed that the trend 
described by Moore’s law will come to an end, and probably 
soon. In contrast, the brain is extremely energy-efficient. With 
many capabilities that are still far beyond modern computers, 
the power of an adult brain is only about 20 watts. Therefore, to 
learn from the brain how to be “greener” is a major goal of BIC. 
With the knowledge obtained in neuroscience, we now know 
that the secret of the brain’s energy efficiency involves various 
factors, including the co-localization of data processing and 
storage, highly distributed processing, and sparse activity. 
Neuromorphic computing aims to implement these features in 
microprocessors, with electronic elements mimicking the 
activities of individual neurons and millions of artificial neurons 
interacting with each other to process information (Merolla et al., 
2014). In the most recent advance in this direction, IBM 
reported that they achieved satisfactory performance in 
complex pattern recognition tasks with a neuromorphic chip. 
Compared with conventional chips, the system reduced the 
energy consumption rate by many orders of magnitude (Esser 
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et al., 2016). It is reasonable to expect that the knowledge 
learned from the brain will enable us to eventually combine 
super computing power with extremely low energy demand in 
the not-so-far-away future. 

The second aspect that the brain can teach us is how to 
achieve better performance in so called cognitive tasks. 
Conventional computers, no matter how powerful, know nothing 
beyond what has been written by their programmers. In addition, 
although they are superfast in crunching large datasets, they 
are incapable of solving multiple tasks that a normal person can 
handle with little effort, such as using language, understanding 
a movie, or driving a car in complex environments. The reason 
behind this capability gap lies in the different ways that 
knowledge/rules are learnt and represented in the system. In 
the brain, the huge amount of knowledge learned by our 
countless ancestors during evolution is stored in the genome 
and expressed in the neural network structure during 
development. Later, through these  well-tuned, highly 
sophisticated networks, more knowledge is gained through an 
individual's interaction with the environment, which is 
represented by hundreds of billions of synapses in the brain 
(Nikolić, 2015). In this sense, compared with the hand-coded 
programs that modern computers rely upon, the brain has a 
much greater capacity to learn and utilize complex rules 
(Baum, 2003). Equipped with a design similar to that of the 
brain and trained by algorithms that allow for highly distributed 
knowledge representation, also like that of the brain, deep 
neural networks—artificial neural networks with many layers 
of processing—have turned out to be very powerful in a 
variety of cognitive tasks, ranging from practical image and 
speech recognition to difficult game play (LeCun et al., 2015; 
Silver et al., 2016). The enthusiasm evoked by such exciting 
advances is enormous across the academic community, 
industry and even the general population. With more 
interactions between neuroscience and machine learning, we 
can be optimistic that the distance from general artificial 
intelligence, at the human level or even beyond it, will become 
increasingly shorter. 

The brain has been the object of modern neuroscience 
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research for more than a century, and artificial neuron networks 
as a tool for information processing were suggested as early as 
the 1940s. So, why is BIC attracting so much attention now? 
On the one side, experimental brain research is at the edge of 
revealing the core principle of the brain. Powerful techniques to 
monitor and manipulate neuronal activities are being rapidly 
applied to both human subjects (noninvasively) and novel animal 
models, including various nonhuman primates and genetically- 
modified organisms. These studies have begun to uncover 1) the 
detailed architecture of brain networks and circuits (e.g., Fan et 
al., 2016), 2) the dynamic rules governing network operation (e.g., 
Yu et al., 2013), and 3) how network and circuit activities give rise 
to motion, perception, and cognition (e.g., Janak & Tye, 2015). 
Such studies provide a solid foundation for BIC. On the other side, 
more efficient algorithms to train artificial neural networks have 
been strengthened by powerful computers, making large, 
complex networks useful for practical purposes (Hinton & 
Salakhutdinov, 2006; LeCun et al., 2015). Thus, exciting 
advances in neuroscience and machine learning, as well as 
rapid improvement in computing power and availability of “big 
data”, have emerged almost at the same time, increasing the 
appeal and value of BIC like never before. Progress in these 
individual areas or in their synergization will no doubt be the 
perpetual driving force behind BIC. 

The design of organisms has provided inspiration for many 
ingenious and elegant solutions in engineering. Now is the time 
to turn our eyes to the pinnacle of biological evolution—the 
brain. Today, BIC is the new challenge for bionics and, in many 
ways, probably the ultimate challenge. 
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ABSTRACT 


Lampreys belong to the superclass Cyclostomata 
and represent the most ancient group of vertebrates. 
Existing for over 360 million years, they are known 
as living fossils due to their many evolutionally 
conserved features. They are not only a keystone 
species for studying the origin and evolution of 
vertebrates, but also one of the best models for 
researching vertebrate embryonic development and 
organ differentiation. From the perspective of genetic 
information, the lamprey genome remains primitive 
compared with that of other higher vertebrates, and 
possesses abundant functional genes. Through 
scientific and technological progress, scientists have 
conducted in-depth studies on the nervous, 
endocrine, and immune systems of lampreys. Such 
research has significance for understanding and 
revealing the origin and evolution of vertebrates, and 
could contribute to a greater understanding of 
human diseases and treatments. This review 
presents the current progress and significance of 
lamprey research. 


Keywords: Lamprey; Nerve; Endocrine; Immune; 
Functional gene 


INTRODUCTION 


Lampreys have seven separated gill holes aligned behind the 
eyes on both sides of the head and appear like eight eyes. 
Lampreys pertain to Cyclostomes, rather than Osteichthyes 
(Kuratani et al., 2002), and are the most primitive agnathans 
(jawless fish) among marine animals, belonging to 
Cyclostomata, Vertebrata, Chordata. Vertebrates are divided 
into the evolutionarily more advanced gnathostomes and the 
inferior agnathans, in which Cyclostomes predominate (Figure 
1) (Janvier, 2006). The existing Cyclostomata contain two 
categories: Petromyzonida and Myxini. The Petromyzonida 
contain Petromyzon in North America, Mordacia in the Southern 
Hemisphere, and Lampetra in China. Lampetra include: 
Lampetra japonica, Lampetra morii and Lampetra reissneri, 
among which Petromyrzon marinus spread to the Great Lakes 
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region, and Lampetra japonica and Lampetra morii to China 
(Freamat & Sower, 2013; Shimeld & Donoghue, 2012). 

Cyclostomata first appeared in the Ordovician, peaked during 
the Silurian and Devonian, and then experienced gradual 
extinction in the late Devonian. The piecemeal discovery of 
fossils has proven that lampreys existed 360 million years ago, 
before the time of dinosaurs, and their shape has stayed almost 
unchanged over hundreds of millions of years of evolution. 
Lampreys are, therefore, universally termed as “living fossils” due 
to their strong resemblance to early fossil material, which reflects 
their important position in vertebrate evolution (Gess et al., 2006). 

Lampreys live in rivers and oceans, and most are migratory 
parasitic. They repeat the cycle of reproduction and death like 
other anadromous fish. The larvae (ammocoetes) spend the 
first 7-9 years of their life in freshwater streams, after which they 
undergo metamorphosis and develop suctorial discs, eyes, and 
dorsal fins, and then migrate to the ocean. When lampreys 
reach full adult size via feeding on fish blood, they become 
sexually mature and migrate back to the streams and rivers 
(Youson & Sower, 2001). 

Rats, mice, Xenopus, zebra fish, globefish, Caenorhabditis 
elegans and Drosophila are well-established model organisms, 
with extensive genome sequencing and scientific studies on 
such species greatly contributing to the health and treatment of 
human diseases (Fuller & Tomé, 2005). Considering their 
special status in vertebrate evolution, studies on lamprey 
species have attracted the growing attention of biologists 
worldwide. With the development of modern molecular 
techniques, lampreys have become one of the most important 
research models for understanding vertebrates. The following 
review discusses current research on the nervous, endocrine, 
and immune systems of lampreys. 
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Figure 1 Phylogeny of select Chordate taxa 
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Labels on the right indicate the names of the monophyletic groupings shown on the left. 


PRIMITIVE FEATURES AND ORGANS OF LAMPREYS 


In marine organisms, lampreys are the most primitive extant 
vertebrates, with the earliest fossil record tracing back to the 
Ordovician. Lamprey research began in the mid-19th century. At 
that time, lampreys were considered as direct ancestors of 
modern vertebrates. Many initial studies focused on comparing 
the morphological differences between adult lampreys and 
higher vertebrates to ascertain their taxonomic positions. The 
lampreys used for those studies were all obtained from the wild, 
with ammocoetes were first captured in 1900 (Reese, 1900). 
With the establishment of in vitro fertilization techniques and 
short-term feeding, it became possible to study lamprey 
embryonic development. Embryonic development of sea 
lamprey was first described by George Piavis (Piavis, 1961), 
with his developmental staging system still widely used today 
(Nikitina et al., 2009). 

The appearance characteristics of lampreys are very 
primitive, and their tissues and organs can represent the initial 
state of vertebrate evolution and development. They do not 
have jaws for predation, and possess median rather than paired 
fins, which are sustained by unsegmented ray-finned cartilages. 


Furthermore, they reserve notochords to support their body axis. 


Their skulls are incomplete, consisting of only a joint of ear 
cartilages between both sides of the ear cartilage capsules, and 
a fibrous tissue membrane over the brain, similar to the skulls of 
higher vertebrates at the early stage of embryonic development. 
The gill cage is composed of gill arches and cartilages, and is 
responsible for supporting respiration. The terminal of the gill 
cage is sealed with pericardial cartilages to protect the heart. 
The bodies and tail muscles of lampreys are also the most 
primitive bow sarcomere of all vertebrates. Their digestive 
systems have not differentiated into stomachs. The blind 
ditches (the original morphology of intestinal villi) are 
differentiated in the intestinal tubes, and the real liver appears. 
The heart, which consists of the sinus, one atrium, and one 
ventricle, appeared in the circulation system of lampreys The 
circulation system is a closed tube type like that of higher 
vertebrates, whereby blood flows in closed vessels. Each part 
of the lamprey brain is arranged in the same plane, without any 
brain partitioning. The mesencephalon does not form a 
bigeminum, and the cerebellum and medulla are not separated. 
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The optic nerve does not develop into an optic chiasma, and 
there are only one or two semicircular canals in the inner ear 
balancer (Janvier, 2008; Kuratani, 2005). 

Considering the features mentioned above, lampreys show 
extraordinarily high primitiveness, indicating that they are 
situated at the lowest level of vertebrate evolution. Research on 
comparative anatomy and embryology reveal that all chordates 
have three main features in common at a certain stage of their 
ontogeny or whole life, that is, the dorsal tube, notochord, and 
pharyngeal gill slits. This phenomenon indicates that all 
chordates originated from a common ancestor As a 
consequence, lampreys are not only pivotal species for 
studying the origin and evolution of vertebrates, but also the 
best model for investigating vertebrate embryonic development 
and organ differentiation. 


LOCOMOTOR BEHAVIOR AND NEURONAL NETWORK 
OF LAMPREYS 


Nerve cells play an important role in an organism's behavior. 
Although studies on various ion channels at the single-cell 
level have achieved rapid development, exploration of the 
relationship between disparate locomotor behaviors and the 
entire neural network operation still face many difficulties. 
Drawing a complex circuit diagram of a neural network 
remains challenging, and the monolithic operation mode of 
the nervous system is hard to decipher with current 
technology. Therefore, the best way to explore such issues is 
to use vertebrates containing simple neural networks. The 
number of lamprey nerves is relatively small, with just one 
brain and one spinal cord, which makes studying their motion 
relatively easy. Hence, lampreys are widely used for studying 
the mechanism of swimming controlled by neural networks. 
The brephic brain structures of different vertebrates show high 
similarity, whereas the adult stage structures exhibit great 
differences. Thus, choosing a brephic brain from a 
representative animal as experimental material can effectively 
clarify the processes of phylogeny and ontogeny of all 
vertebrate brains. In the evolution of chordates, the basic 
structure of the brain has not appeared in amphioxi, a close 
relative of vertebrates. Instead, amphioxi have cerebral 
vesicles, which are slightly dilated lumens at the front of the 
neural tubes (Holland & Holland, 2001). The brains of 


lampreys already possess the basic structures of the 
vertebrate brain, which are divided into the forebrain, 
diencephalon, midbrain, cerebellum, and medulla (Murakami & 
Kuratani, 2008). In embryonic development, biological 
ontogeny will repeat its phylogeny, which means that the 
development process of animal embryos usually repeats the 
important stages of germline evolution. The ontogeny of 
lamprey brains indeed replays the phylogeny of vertebrate 
brains. The structure of the lamprey brain can reflect the 
elaborate structural changes of higher vertebrate brains. Thus, 
as lower vertebrates, lampreys are one of the most valuable 
animal models for the study of vertebrate brains. 

Lampreys have many visible neurons, making them highly 
suitable for microelectrode experiments. Moreover, their 
neurons can survive for a long time in-vitro in separation 
devices, which is helpful for bionics (Murakami & Watanabe, 
2009). Professor Grillner and colleagues (Biró et al., 2008; 
Ménard & Grillner, 2008; Saitoh et al., 2007) conducted 
extensive in-depth work on the mechanism of neural network 
control of lamprey movements. They described the specific 
structure of the nervous system, and provided a theoretical 
description of the activity patterns of neurons in the network 
structure. They emulated the swimming motion 
characteristics of lampreys by utilizing neuronal activity 
described with the Hodgkin-Huxley (HH) model, which 
contained ion channels and included 3D models to simulate 
forward and backward swimming motion and changes in 
swimming direction. Kamali et al. (2013) constructed an 
experiment in which a simulated lamprey swam around a 
virtual environment containing visually detectable objects, 
which addressed mechanisms for the generation and 
selection of visual behaviors in anamniotes. Ericsson et al. 
regulated dopamine levels to study how the strength of 
synaptic connections affected mode movements in lampreys 
(Ericsson et al., 2013). Researchers have also investigated 
the effects on lamprey movement rhythm when their 
cephalic neuromeres and spinal neurons suffer from 
external stimuli (Cinelli et al., 2014; Kinkead, 2009). 

Rétaux et al. (Rétaux & Kano, 2010) studied Sonic 
Hedgehog/Hedgehog (Shh/Hh) signaling in lampreys, and 
discovered that embryonic midline signaling is the driving force 
of vertebrate forebrain evolution. As the most primitive 
vertebrate, the neural crests and structures of lampreys are 
differentiated from the neural crests of jawed vertebrates, 
although they still lack some important neural crest derivatives, 
such as sympathetic ganglions. Research further suggests that 
the migration path of lamprey neural crest cells is similar to that 
of higher vertebrates. 

These studies again show that lampreys are excellent 
animal models for studying the evolution of the vertebrate 
nervous system. The nervous system, neural synapses, and 
synaptic neurotransmitters of lampreys represent the most 
primitive state in vertebrates. Uncovering the mysteries of 
the human brain is a common goal of neurobiologists 
worldwide. Lampreys, as the headwaters of vertebrate 
evolution, might hold the key to unlocking human cognition 
and the brain maze map. 


LAMPREY ENDOCRINE SYSTEM 


Over the past decade, molecular genetics have developed 
rapidly. Gonadotropin releasing hormone (GnRH) is a key 
signaling molecule of the hypothalamic-pituitary-gonadal axis. 
In lampreys, olfactory stimuli may play a major role in regulating 
GnRH secretion. Research shows that GnRH in lampreys is 
divided into GnRH-l, GnRH-Il, and GnRH-IIl, and is mainly 
distributed in the hypothalamus, where it controls the pituitary 
gland and thus plays an important role in gonadal development 
and maturation of gametes (Decatur et al., 2013). The lamprey 
endocrine system also plays a part in metamorphosis. In the 
initial metamorphosis stage of parasitic and non-parasitic 
lampreys, the levels of thyroid hormone and liothyronine in 
serum drop sharply (Freamat & Sower, 2010). In fact, the 
compounds that can affect thyroid hormone synthesis, such as 
KCIO, and goitrogens, are capable of initiating metamorphosis 
in lampreys, and also lowering the level of serum thyroid 
simultaneously (Manzon & Youson, 1999). Sexual maturity in 
lampreys is slow, which does impact its use as a model 
organism. It is not only beneficial for lamprey research, but also 
helpful in finding therapies for human diseases caused by 
hormonal disorders if manual methods to accelerate the 
lamprey maturation process are used. 

The hypothalamus of higher vertebrates produces corticotrophin- 
releasing hormone (CRH) under stress. This hormone can 
stimulate the anterior pituitary cells to release adrenocorticotropic 
hormone (ACTH), with adrenal glands then induced to release 
cortisol or corticosterone. However, there are no such steroid 
substances in lampreys. Roberts et al. (2014) confirmed that 
there was one potential corticosteroid, 17, 21-dihydroxypregn-4- 
ene, 3, 20-one, that participates in the endocrine regulation of 
compression resistance in lampreys. The mode of endocrine 
changes in lampreys under stress might be the prototype of 
how the endocrine system varies in higher vertebrates when 
they face ambient stresses. 

It is not only advantageous for scientific research to find a 
way in which to promote the growth of lampreys, but studies on 
their endocrine systems can also help capture the lampreys 
which usually cause fishery damages. Due to their parasitic 
habits, lampreys can have considerable impact on fisheries. 
Fish farmers will disrupt lamprey breeding, development, and 
migration in order to prevent damage to their fisheries. For 
example, the use of the drug 3-trifluoromethyl-4-nitrophenol 
(TFM) can kill lamprey larvae, while avoiding other fish and 
organisms (Birceanu et al., 2014). Because of the high cost of 
TFM, however, it is imperative that new drugs are developed to 
mitigate lamprey impact. Lamprey hormones might be a 
suitable choice, whereby mature lampreys can be prevented 
from breeding by using specific hormones or derivatives to lure 
them to a predetermined fishing spot. 

The lamprey endocrine system is simpler than that of humans, 
and the endocrine organs and related molecules are more 
primordial. Comparing the histology and biochemistry of the 
primitive lamprey endocrine system with the well-differentiated 
human endocrine system has enabled researchers to better 
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comprehend the evolution of the vertebrate endocrine system. 
LAMPREY IMMUNE SYSTEM 


The origin and evolution of the adaptive immune system (AIS) 
remains an important field in immunological research. Jawless 
vertebrates are generally considered to have evolved the AIS. 
Therefore, clarifying the adaptive immunity of Agnatha is of great 
significance in revealing the origin and evolution of the AIS. Early 
experiments showed that lampreys could produce specific 
agglutinin memorially contraposing particulate antigens, such as 
anthrax and sheep blood cells (Fujii et al., 1979), indicating that 
lampreys should have an AIS. Initially, however, T-cell receptors 
(TCR), B-cell receptors (BCR), and major histocompatibility 
complexes (MHC), which are symbolic molecules of the AIS, 
could not be tested in lampreys and was not clarified until an 
alternative AIS set was found (Cooper & Alder, 2006). 

Adaptive immunity relies on the formation of lymphocytes. 
Scientists have found lymphocyte-like cells in lampreys that 
are similar to the lymphocytes in higher vertebrates. Shintani 
et al. (2000) found homologous genes of PU.1/Spi-B (purine 
box 1/spleen focus-forming virus integration B), a gene 
involved in the differentiation of lymphocytes in jawed 
vertebrates, expressed in the lymphocyte-like cells of the 


digestive tract and inexplicably the ovary in lampreys. Thus, 
these results validated the existence of lymphocytes in the 
jawless vertebrates. 

To find the molecules that mediate the immunoreaction of 
jawless animals, Pancer et al. (2004) identified a large number 
of cDNA clones encoding multiple leucine-rich repeat (LRR) 
modules, and instead of the V (D) JC composition of 
immunoglobulin-based BCRs and TCRs observed in jawed 
vertebrates, variable lymphocyte receptors (VLRs) of jawless 
vertebrates were found. VLRs are composed of variable LRR 
sequences, and are produced by novel molecular 
rearrangement. The germline VLR gene does not have an 
intact structure to encode any proteins. Many LRR-encoding 
cassettes with highly diverse sequences are adjacent to this 
gene. During the development of lymphocyte-like cells, these 
LRR-encoding cassettes are chosen randomly and sequentially 
incorporated into the germline VLR gene. Not only is the 
sequence of each LRR cassette diversified, but the number of 
cassettes inserted into the VLR gene is variable (1-9). Thus, a 
computational assessment of mature VLR gene sequences and 
LRR cassettes used to assemble them can predict a potential 
repertoire, which is almost equivalent to that of Igs (>10"). 
Such diversity is enough to identify any antigen outside (Hirano 
et al., 2011) (Figure 2). 
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Figure 2 Rearrangement of antigen receptors of lampreys 
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In the vicinity of the incomplete germline VLR gene are multiple LRR-encoding cassettes. During the development of lymphocyte-like cells, these 


cassettes are sequentially incorporated into the VLR gene. LRR1, the first 
and LRRVe, the terminal 24-residue LRR with a distinct sequence signature. 
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18-residue LRRs; LRRV, 24-residue LRRs, whose number is variable; 


Monoclonal antibodies have been widely used in research 
and clinical practices. However, immune globulins are costly, 
the manufacturing techniques are intricate, and the 
physicochemical properties are unstable. These defects have 
prompted research on finding alternatives. The discovery of 
VLRs affords a new perspective for the development of 
antibody reagents or drugs. The molecular weights of VLRs are 
about 15-45x10?, while the relative molecular mass of IgGs can 
reach 156x107. Small molecular weights are more suitable for 
in-vivo transport of drugs and penetration into blood vessels, 
which means that VLRs are more applicable for the 
development of antibody drugs (Boehm et al., 2012). The 
monoclonal VLRs prepared by Herrin et al. (2008) had 1 000- 
fold higher affinity than IgGs at the same concentration, and 
could only be eluted under alkaline conditions (pH>11). In 
addition, the relationship between lampreys and humans is not 
close, so VLRs have good specificity as diagnostic reagents, 
eliminating cross-reactivities. The molecular structure of VLRs 
is relatively simple, with only one peptide chain. Compared with 
the hybridoma techniques of preparing monoclonal antibodies, 
the manufacturing techniques of VLRs are easier and 
production costs are lower, which makes it convenient for large- 
scale industrial production. 

The discovery of immune organs is still a hot topic in lamprey 
research. As confirmed recently, lampreys do not have a 
thymus, lymph nodes, spleens, or tonsils, nor do they develop 
any specific immune organs and tissues. However, they must 
possess a powerful immune function to have withstood various 
environments and the invasion of germs over hundreds of 
millions of years. Accordingly, lampreys must contain a powerful 
immune organ to support the immune function. In lamprey 
larvae, researchers have found that lymphocytes might arise in 
the typhlosoles (i.e., sunken area of the intestinal epithelium) 
and renal folds (Bajoghli et al., 2011). However, information on 
the exact immune organs of adult lampreys remains uncertain. 
Recent lamprey studies by our group found a strip of medullary 
tissue throughout their entire trunk and tail above the nerviduct 
(Pang et al., 2015). It was preliminarily confirmed that this tissue 
was mainly made up of fat cells, hematopoietic stem cells, and 
lymphocyte-like cells. We termed it the “immune body" due to 
its potential immune function (Figure 3). In the cDNA library of 
the immune body constructed by our group, we found abundant 
immune molecules that still require further study. It is possible 
that a very large immune organ, together with substantial 
immune molecules, has allowed lampreys to survive till the 
present day during despite the impact of environmental 
migration and bacterial infection during evolution. 


immune body 






oral gland gills cage 


Figure 3 Location of immune body, oral gland and gill cage of lamprey 


It has been a widely held belief that adaptive immunity in 
higher vertebrates evolved from the innate immune system (Du 


Pasquier, 2005). However, the discovery of VLR and its 
restructuring mechanism in jawless vertebrates impacted the 
innate immune evolutionary theory. This type of LRR-insert 
restructuring mechanism is likely the evolutionary origin of 
adaptive immunity. Currently, human immune-related diseases 
have become the culprits of health injury, even death. The in- 
depth exploration of the lamprey immune system is necessary 
for defeating serious human illnesses such as cancer, AIDS, 
and autoimmune diseases. 


GENOME AND FUNCTIONAL GENES OF LAMPREYS 


Since lampreys are living fossils of ancient animals, the lamprey 
genome could be described as a molecular fossil of ancient 
genetic information. Therefore, study of the lamprey genome 
could help clarify the origin and evolution of vertebrates. The 
lamprey genome is a record of the genetic information of 
vertebrates hundreds of millions of years ago, and exhibits 
many novel genes and proteins. Potter et al. (1968) discovered 
that chromosome number showed great variation between 
different sub-species of lamprey, 2n= 60-168. Because the 
number of lamprey chromosomes is large and the volume of 
each chromosome is very small, karyotype analysis can be 
difficult. However, it does provide some evidence for the 
chromosome doubling theory of vertebrates. The above 
phenomenon illustrates that the lamprey genome still remains 
in an intermediate stage of evolution from invertebrates to 
vertebrates. 

Over the years, studies on lampreys have been confined to 
ecology and physiology, while genome and proteomic research 
remains scarce, especially for some important functional genes. 
Our group has conducted extensive lamprey genomic and 
proteomic research. We successively constructed lamprey 
cDNA libraries of salivary glands, liver, blood and immune body 
and other tissues, and sought out many novel genes and 
proteins with important physiological functions. For example, we 
found a Bruton's tyrosine kinase like (BTK) gene expressed in 
lamprey immune-related tissues, which might be involved in the 
immune response of lampreys (Wu et al., 2012). Reactive 
oxygen species modulator | (Romo |) and Peroxiredoxin (Prx) 
homolog were discovered, and are capable of inducing the 
production of reactive oxygen, resisting peroxidation of 
erythrocytes, and adjusting the cellular oxidation-deoxidation 
environment (Sun et al., 2010; Zhao et al., 2013). Transcription 
factor NFKB inhibitor I-KB-£ was also found in lampreys, and its 
expression was reduced under the stimulation of LPS (Su et al., 
2013). Among the many functional genes found, the genes and 
proteins in the oral gland are special. In the evolutionary 
process, lamprey secreta could prevent blood coagulation of 
host wounds. Lampreys have a special way of life, which results 
in many proteins and genes with distinctive functions. The oral 
gland is wrapped in a special colored sac in the hypobranchial 
musculus (Figure 3), and can secrete an eel plasmin substance 
to prevent blood from clotting. We also found other functional 
proteins in lamprey oral glands, including those containing RGD 
(Arg-Gly-Asp) motif disintegrin-like proteins, retinoic acid, and 
interferon-induced lethality protein-19 (GRINM-19), tumor 


Zoological Research 37(5): 263-269, 2016 267 


translationally controlled protein (TCTP), and gelsolin-like 
protein (Sun et al., 2008; Wang et al., 2010; Xiao et al., 2012). In 
subsequent future studies, we are hopeful to exploit intravenous, 
oral, or gene therapy drugs used for anti-tumor, anti-thrombotic, 
anti-virus, hypoglycemic, analgesic, and  anti-heart failure 
treatment on the basis of these special function proteins. 

Studies on lamprey functional genes can help to understand 
the detailed biological evolutionary process of different genes, 
and can also contribute to finding associated proteins or peptide 
drugs, which might be favorable for human disease therapy. 
Therefore, whether from the perspective of scientific research 
or application, studying the function of lamprey genes has far- 
reaching significance. 


CONCLUSIONS AND PERSPECTIVES 


Due to their unique way of life and appearance, lamprey 
research began as early as the 19" century. Many studies 
focused on comparative biology, especially in the field of early 
vertebrate evolution. Lampreys have also been used in the 
study of evolutionary embryology and developmental biology. 
With the emergence of modern molecular techniques, lampreys 
have become one of the most important models for the 
exploration of vertebrate origins. 

Ancient lamprey fossils have provided evidence of their 
emergence more than 360 million years ago. Their external 
morphology shows no obvious differences compared with 
existing lampreys, which is known as "evolutionary diapause". 
Lampreys exhibit primitive external morphology and physiology 
structure: no paired appendages; only horny teeth formed by 
the epidermis; notochord retained throughout life; skull 
incomplete, without top, equivalent to the early stages of 
embryonic development of the vertebrate skull; less muscle 
differentiation, original arrangement of sarcomeres; lower 
degree of brain development layer, and two semicircular canals. 
These original characteristics suggest that lampreys are the 
evolutionary node from invertebrates to vertebrates. Thus, due 
to their critical species status, lamprey research is fundamental 
to biological evolution and genomic studies. 

The rapid development of cell biology and biochemistry 
technology, especially genomics and proteomics, has pushed 
lamprey research to new heights. Lampreys have a powerful 
and unique immune system and many valuable genes. Based 
on lamprey genome research, important human disease-related 
genes might be discovered. In addition, the pathogenesis of 
human diseases can be revealed by establishing lamprey 
genetic models of human diseases, which could provide the 
theoretical basis for in-depth study of the pathogenesis and 
treatment of diseases and development of new drugs. In recent 
years, lamprey numbers have drastically reduced, so lamprey 
research is important. 

As a model animal for the evolution and development of 
vertebrates, lampreys have high research value due to their 
connection to vertebrates and invertebrates in evolution. 
Lampreys are not only key species for the exploration of the 
origin and evolution of vertebrates, but also the best model for 
studying vertebrate embryogenesis and organ differentiation. 
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The construction of lamprey models will occupy a central 
position in vertebrate evolution and developmental studies 
during this century. 
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ABSTRACT 


In recent years, wildlife conservation has attracted 
great public attention. However, substantial 
distinctions can be found in the prevailing concepts 
of wildlife conservation, particularly with the recent 
notion that emphasizes animal rights. Wildlife 
welfare and wildlife rights are not synonymous, with 
welfare more compatible with the reasonable and 
legal utilization of wildlife. The key to scientific 
wildlife conservation is the appropriate awareness 
and appreciation of the relationship between wildlife 
conservation and utilization and the theoretical basis 
of holism. Nevertheless, rational biases regarding 
the public’s understanding of wildlife conservation 
and the spread of information via social media still 
exist. As such, expansion of the concept of scientific 
wildlife conservation requires the application of 
several measures. Wildlife conservation researchers 
should be regarded as the most important 
disseminators of scientifically-based information, 
with education in schools and universities of growing 
importance. Furthermore, the media should shoulder 
the social responsibility for the accurate 
dissemination of conservation information. 


Keywords: Wildlife conservation; scientific concept; 
Media; Dissemination. 


INTRODUCTION 


Wildlife conservation has two meanings. One is the 
preservation of both species and species diversity, the other is 
based on animal welfare, which is primarily aimed at wildlife in 
captivity (Lu, 2009). Conservation education is an important 
component of environmental education, and is aimed at 
expanding human awareness of conservation biodiversity and 
at changing environmental attitudes and behaviors to promote 
conservation through education and practical activities. 

Wildlife conservation education forms part of conservation 
education. Since environmental concerns have increased 
across all society, wildlife conservation has become a 
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significant social issue. However, there are considerable 
differences in the concepts of wildlife conservation, with several 
plausible protection ideas currently debated. Some people 
believe that wildlife conservation should incorporate the 
protection of all animals, and follow the principles of animal. In 
addition, such protection includes resistance to live bear bile 
used in traditional Chinese medicine and prohibition of hunting 
and related activities. Furthermore, absolute conservation has 
strong public sensibilities, which can result in extreme wildlife 
conservation activities, thus welfare and animal rights. In 
contrast, others believe that wildlife conservation should be 
based upon scientific attitudes and strategies. Unfortunately, 
absolute protection currently dominates public opinion and 
sympathy (Zhang et al, 2015), with scientific discussion and 
rational thinking often discounted. ' 

Absolute protection includes the random release of animals, 
which has led to the invasion of alien species weakening and 
hindering the process of wildlife conservation itself (Karanth et 
al., 2008). 

Since public attitude plays a very important role in the 
formulation and implementation of wildlife conservation 
management policies, some of the main problems that need to 
be solved in the construction of an ecologically-based 
civilization are defining wildlife conservation, promoting 
scientific wildlife protection, and encouraging an objective view 
of wildlife conservation by the public. 


CONCEPT OF SCIENTIFIC WILDLIFE CONSERVATION 


To understand the public’s attitude to wildlife conservation, we 
must first understand the relationship between humans and wildlife. 


Relationship between humans and wildlife 

Wildlife have been an important natural resource during the 
process of human evolution, and animal farming has supplied 
the necessary material for long-term human development. Thus, 
the relationship between wildlife and humans is close and 
complicated. 
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Wildlife as an important resource in human evolution 
Research has shown that Homo habilis and Homo sapiens 
(archaic) increased the proportion of meat protein in their diets 
by hunting, resulting in the rapid development and increased 
volume of the brain (Wang, 2004). In addition, meat was crucial 
to the reproduction and evolution of prehistoric humans (Niche, 
1995). Compared with Homo habilis, Homo sapiens (archaic) 
used fire to cook and soften food, resulting in large changes in 
human morphology (Qiao, 2011). The shortening of food 
chewing and digestion time resulted in significant changes in 
mandible and dental morphology, such as jaw retraction, 
smaller teeth, rostrum retraction, and smaller crinium viscerale 
(Yu & Zhai, 2004). In addition, human learned to use fire and 
the consumption of meat protein, which is more easily absorbed 
than raw meat, following the introduction of fire and cooking 
promoted the development of the human brain and body. 
Furthermore, hunting activities improved human cognition and 
problem solving abilities, and promoted the evolution of physical 
fitness (Qiao, 2011) In short, the consumption of meat 
promoted human body health by reducing disease and 
strengthening the functions of the brain and other organs (Qiao, 
2011). 


Importance of wildlife in the long-term development of 
human society 

Wildlife have played an important role in the development of 
clothing, medical materials, experimental models, and 
scientific research. Animal skins and fur have been used for 
clothing for millennia, and even today are symbols of fashion. 
Wild animals have also provided nutrition for humans, forming 
a significant proportion of our diet. The domestication and 
farming of wildlife, the advancement of feed technology, and 
the invention of meat and milk production can be considered 
the three revolutions of the human diet. These developments 
improved the fat and protein proportions in the diet, and 
enhanced the development of the human body and brain. At 
the early stage of Homo sapiens, wide-spread farming and 
preliminary formation of agricultural production increased 
human food abundance and variety, which again contributed 
to an increase in brain volume. At the late stage of Homo 
sapiens, with the arrival of primary agriculture, animal 
husbandry, and the industrial revolution, the structure, habit, 
and concept of the human diet has tended towards stability. 
Correspondingly, the physical form has remained predominantly 
unchanged (Qiao, 2011). 

Medical advancement has also been significant in the 
development of humans. Animals form a vital part of traditional 
Chinese medicine, as both raw materials and secretions, and 
play an important role in the continuation of Chinese 
descendants (Liu, 2014). Wildlife also provide raw materials for 
scientific and medical research. 

Humans and wildlife are closely interrelated, especially in 
regards to culture, traditional Chinese medicine, food, hunting, 
and eco-tourism. Therefore, it is not easy and appropriate to 
separate wildlife from human. It is both socially and scientifically 
important, therefore, that we clarify how best to protect wildlife 


and upon which concept of wildlife conservation this protection 
is based. 


Differences between wildlife welfare and rights 

Wildlife rights and welfare are two different aspects of wildlife 
conservation. Those that support wildlife rights argue that 
animals possess thoughts, desires, consciousness, and 
memories, and feel emotion and pain the same as human 
beings. As such, they have the right to not being hurt or 
exploited, and it is our fundamental obligation to not use 
animals in research or merchant farming. Professor Tom Regan, 
an American philosopher in animal rights theory, states that 
animals are also the subject of life and, like humans, possess 
both heart and psychology, and are therefore deserving of 
rights (Cai, 2006). Many animal rights organizations oppose 
animal experiments. For example, in March 2012, China 
Southern Airlines canceled a shipment of crab-eating 
macaques (Macaca fascicularis) from Guangzhou to Los 
Angeles because of PETA's (People for the Ethical Treatment of 
Animals) protest. Some people think that we should be kind to 
animals, and should stop conducting animal-based 
experimental research (Wadman, 2012). Conversely, those that 
support animal welfare are concerned with the rational and 
humane use of wildlife. They state that humans are morally 
superior to wildlife, and thus should be allowed to use animals 
for their own benefit, though without unnecessary pain (Cai, 
2006). Professor Carl Cohen, who argues against animal rights, 
believes that animals do experience pain and therefore must be 
care for, but that wildlife and animal-based experiments should 
be conducted, particularly in regards to the relief of human pain 
and disease (Huang, 2014). 

However, what are the fundamental differences between 
wildlife rights and welfare? At present, more than 100 countries 
and regions have established animal welfare regulations, 
especially within European and American countries where 
detailed and interoperable animal welfare legislation started and 
developed earlier. For example, the mistreatment of animals, 
poor breeding conditions, and abandonment of pets are likely to 
be prosecuted, and the criminal to face legal punishment in 
such countries. Even in legitimate experimental institutions, 
poor and or illegal animal treatment will result in possible legal 
action (Wang, 2009). However, this does not mean that in these 
regions, animals have the same legal protection as human 
beings. Thus, rights and welfare differ in definition and 
application, and it is unrealistic and unscientific to advocate for 
wildlife rights. 


Differences between absolute and scientific wildlife 
conservation 

Absolute wildlife conservation 

Although the differences between wildlife rights and welfare are 
very clear, in reality extreme events always exist in wildlife 
conservation. For example, the random release of the red- 
eared turtle (Trachemys scripta elegans), which led to a 
substantial reduction in China's native turtles (Zhou, 2010), and 
the extreme publicity of vegetarianism. These people protect 
any animal without reason or understanding of the 
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consequences. There are two reasons for this extreme situation. 
The first is exaggerating or confusing wildlife rights with wildlife 
welfare, in which animals and humans are regarded as equals 
and animal use of any kind is considered unethical. The second 
is the imposition of a personal choice on the public. For 
example, some people are vegetarian or regard animal 
companions as family members, and therefore oppose the use 
of animals based on their personal attachment. This concept of 
absolute conservation could misleadingly influence the attitude 
of the general public in regards to wildlife protection. Conversely, 
the concept of scientific conservation is the dialectical and 
objective understanding of the relationship between wildlife 
protection and utilization. Wildlife protection and utilization are 
inseparable, and their division could impact the balance 
between humans and nature. Scientific conservation requires a 
rational view of animal protection and the scientific utilization of 
wild animals (Zhang, 2006). It is not rational to oppose all 
activities in which animals are used in the name of animal 
protection, nor equate wildlife welfare with wildlife rights 
because of a love for animals. As to personal attitudes towards 
animals, it is not appropriate to impose or force a personal 
choice on the public, resulting in distortion of public attitude to 
wildlife conservation. 


Wildlife conservation based on holism 

The key to the concept of scientific wildlife conservation is to 
understand the relationship between the protection and 
utilization of wildlife based on holism, which asserts that a 
system is an organic whole and each part cannot be 
understood separately. From the perspective of holism, the 
protection and utilization of wildlife is unified, not conflicted. 
However, wildlife protection and utilization are often incorrectly 
thought to be contradictory. 

The most typical example is wildlife hunting, in which those 
that oppose it argue that animals still require protection. 
However, the reasons why people disagree with hunting are 
based on a general lack of in-depth understanding. The 
principle of wildlife management states that disordered and 
arbitrary hunting will cause serious damage to wildlife, whereas 
methodical and well-managed hunting activities can be 
beneficial to wildlife conservation (Cao et al, 2014). 

Well-organized hunting not only has ecological benefits, but 
also economic and social ones. From the perspective of 
ecology, ordered wildlife hunting can help regulate a population. 
For example, the poor ecological balance resulting from the 
over-population of herbivores due to low numbers of large 
carnivores can be stabilized through the artificial control of 
herbivore populations. 

Wildlife hunting in many countries has proven that ordered 
hunting can have a positive effect on maintaining the ecological 
balance (Cao et al, 2014), as well as economic and social 
benefits. According to the United States Fish and Wildlife 
(USFWS) Service, the effect of wildlife hunting on economics 
was over $65 billion in 2002, and direct income from wildlife 
hunting was $22 billion (Cao et al, 2014). Hunting also offers 
more than 700 employment opportunities for people, with one in 
every 15 Americans engaged in hunting activities. If these 
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individuals are considered as a company, their tax ranks as 35 
in the USA fortune 500 companies. Furthermore, the 2011 
USFWS report states that 1.37 million people over 16 years 
participated in hunting activities, with total hunting expenses of 
$33.7 billion (Tang & Zhou, 2013). 

Another example of typical absolute conservation is the 
random release of wildlife. To protect animals, some 
organizations and individuals have planned and instigated the 
random release of wildlife into the field, even buying such 
wildlife for liberation. This behavior not only encourages the 
illegal trade of wildlife, but can also lead to the destruction of the 
natural environment due to the invasion of foreign species 
(Zhou, 2010). Although such negative results often deviate from 
the original intention, they do highlight the need for professional, 
rather than fervent and inexpert, protection. 

Scientific and professional wildlife conservation is, therefore, 
very different from extreme conservation. Protection and 
utilization of wildlife are not separate entities. The relationship 
between wildlife conservation and exploitation needs to be 
treated dialectically. The opposite of protection is the 
destruction of wildlife and the resulting ecological imbalance, 
not the limited use of certain wildlife products. 


INADEQUACIES IN THE MEDIA AND PUBLIC REGARDING 
WILDLIFE CONSERVATION 


With the expansion of wildlife conservation publicity after the 
passing of Wildlife Protection Law in 1988 (http://www. 
china.com.cn/chinese/law/647629.htm), wildlife conservation 
has gained significant public attention. However, due to the lack 
of scientific and objective understanding of wildlife conservation, 
social media are more focusing on how lovely the wildlifes are, 
the urgency of protection and the destructiveness of overusing. 
Social media in China does not systematically introduce 
scientific strategies of wildlife conservation to the public, which 
helps explain why people easily hold to the concept of absolute 
wildlife conservation. Slogans such as “stop the use, will stop 
the killing", “plead for animals", “animals and humans are equal”, 
and "we are family’ have gone beyond the fundamental 
purpose of wildlife protection (Zhang, 2013) and have pushed 
the relationship between humans and wildlife to the extreme. 

The food chain tells us that the killing of animals is usual in the 
wild. Without killing, neither carnivores nor herbivores can survive. 
Eating or being eaten forms part of the ecological balance. From 
this perspective, wildlife conservation is the protection of the 
species, rather than the protection of the individual. 

If absolute wildlife conservation is non-scientific and non- 
rational, the direction and scientific process of wildlife 
conservation can be hindered, and the balanced relationship 
between humans and wild animals destroyed, resulting in 
ideological confusion in the field of wildlife conservation. 


DISSEMINATION OF SCIENTIFIC WILDLIFE CONSERVATION 
Popularization of scientific wildlife conservation by 


researchers 
It is difficult for many people to obtain detailed information on 


environmental issues, often relying on social media, which can 
be biased and result in public misunderstanding (Coyle, 2005; 
Ladle et al., 2005). Scientific researchers are responsible for 
encouraging the public to understand the concept of and 
requirements for biological diversity conservation (Sunderland 
et al., 2009). As such, communication with the public must be 
effective, which requires active input from wildlife researchers 
and the development of effective methods to communicate with 
different people in regards to conservation education. 

Sometime, people are more willing to believe in something 
that they think increases their personal values, whereas, ignore 
scienfic data. One possible solution is to help people to find out 
the root of problem and establish their values of environmental 
protection on the basis of scientific understanding (St Clair, 
2003). An excellent scientific communicator understands that 
scientific data are often not effective at persuading the public, 
with more specific and targeted strategies required. 
Conservation biologists should integrate the local culture into 
the overall planning of scientific communication, and not try to 
change behavior and thinking by data alone. 


Conservation education in schools 

Research has shown there are connections between high-level 
knowledge and high-level action for environmental protection 
(Rickinson, 2010). Although environmental education has 
demonstrated a growing trend in certain developed countries, 
the education is still limited due to poor access to educational 
courses (Cutter-MacKenzie & Smith, 2003). 

Li et al. (2014) investigated the concept of conservation by 
comparing changes in the attitudes of university students 
towards wildlife after a professional lecture about sustainable 
conservation. According to the research, although the grade 
and hometown of students had certain influences on their 
attitudes to wildlife conservation, the expert lecture had a 
positive effect on students understanding (Li et al., 2014). Thus, 
attitudes towards wildlife conservation are variable, and not 
necessarily scientifically based. Research has stated that 
wildlife conservation education for university students should 
include the concept of scientific wildlife conservation through a 
variety of practical platforms and curricula (Zhang et al., 2015), 
which can be incorporated into existing courses in schools and 
universities and spread through an interdisciplinary approach 
(Sodhi et al., 2003). Publicity and education in schools can 
effectively enlarge the range of people interested and educated 
in the scientific protection of wildlife. It is not only the 
responsibility of universities to offer sustainable wildlife 
conservation courses, but also that of elementary and junior- 
senior high schools, which can lay a scientific foundation for 
more comprehensive education and sustainable wildlife 
conservation. 


Popularization of scientific wildlife conservation by the 
media 

Formal educational institutions are not the only channel from 
which the public can obtain their environmental knowledge. 
Social media and relationships are also important. In particular, 
social media can be considered as a "school of the air" where 


people can gain knowledge on environmental protection (Ors, 
2012). Generally speaking, social media is more capable than 
researchers in delivering information. However, the frequent 
misunderstanding of scientific knowledge by the media can lead 
to unsatisfactory or inaccurate reports (Pace et al., 2010). This 
is especially evident in the field of wildlife conservation. 
Nowadays, the concept of absolute wildlife conservation 
dominates social media, making rational discussion difficult. 
Because social media plays an important role in information 
distribution, it is very important that the concepts of scientific 
wildlife conservation are clear. In particular, researchers should 
be incentivized and willing to popularize their scientific data, in 
addition to publishing their articles in academic journals, and 
articles published in social media should be emphasized. Thus, 
the public could gain scientific knowledge from media, such as 
newspapers and blogs (Sodhi, 2007), while scientists expand 
access to their research. 

In summary, public awareness and knowledge of the concept 
of scientific wildlife conservation is crucial for animal protection, 
and effective dissemination of scientific wildlife conservation 
has become a key issue in the field. The successful spread of 
information relies on the devotion of time and energy to 
overcome the obstacles caused by different cultures, 
educational backgrounds, beliefs, economic conditions, and 
regions. 
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ABSTRACT 


Biological invasion is a pervasive negative force of 
global change, especially in its effects on sensitive 
freshwater ecosystems. Even protected areas are 
usually not immune. Ptychobarbus chungtienensis is 
a threatened freshwater fish now almost confined to 
Bita Lake, in the Shangri-La region of Yunnan 
province, China. Its existence is threatened by the 
introduction of non-native weatherfishes (Misgurnus 
anguillicaudatus and Paramisgurnus dabryanus) by 
an unusual method known as ‘prayer animal release’. 
Periodic surveys revealed the ratio of invasive 
weatherfishes to P. chungtienensis has been 
increasing since the former species was first 
recorded from the lake in August, 2009. 
Ptychobarbus chungtienensis shows low genetic 
diversity in the relict Lake Bita population. 
Weatherfishes, however, have highly successful 
survival strategies. The degree of dietary overlap 
between the species is alarming and perhaps critical 
if food is found to be a limiting factor. 


Keywords: Biological invasion; Threatened fish; 
Prayer animal release; Genetic diversity; Dietary 


INTRODUCTION 


Biological invasions have become a pervasive agent of global 
change not only in homogenizing the earth's biota (Lodge, 1993; 
Simberloff et al., 2013), but also in challenging the conservation 
of biodiversity and natural resources (TEEB, 2010). Freshwater 
ecosystems are particularly sensitive to such invasions (Villéger 
et al., 2011), where invasive species become one of the major 
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threats to the global freshwater biodiversity (Geist, 2011; 
Magurran, 2009). Although without enough attention, China 
is one of the countries suffering from serious biological 
invasions (Xu et al., 2006; Xing et al., 2016). When taking 
the intentional transplantation (Kang et al., 2013), escape of 
domestication stocks (Yang et al., 2011), and accompanied 
by other forms of introduction (Chen, 2010) into 
consideration, the biological invasion in the freshwater 
ecosystem in China is quite severe. Generally speaking, 
lacustrine ecosystem (lakes, ponds, etc.) are more sensitive 
than riverine ecosystem (rivers, streams, etc.), and the 
extinction, replacement of indigenous fish by introduced 
species is extremely common all over the lakes of China 
whether from plain lakes in the east (Chen et al., 2010), or 
plateau lakes in the west (Chen et al., 1998). 

Bita Lake is a plateau freshwater lake in the Yunnan Bitahai 
Nature Reserve, which is a major part of the Putatso National 
Park. It is also the core area of the Bitahai wetland, one of the 
Wetlands of International Importance (Ramsar sites in 2004), 
with an elevation of about 3 568 m a.s.l. and a catchment area 
159 hm? (Yang & Ji, 2010) (Figure 1A-C). The fish fauna in Bita 
Lake is extremely simple, Ptychobarbus chungtienensis, one of 
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the plateau schizothoracinae fishes (Cyprinidae), was recorded 
as the only fish species (Figure 1D), which was considered to 
be extraordinarily unusual in the evolution of species and the 
formation of lakes (Chen & Chen, 2002). Ptychobarbus 
chungtienensis is a threatened freshwater fish that has been 
listed as Endangered on the China Species Red List (Pan et al., 
2010; Wang & Xie, 2004). It was originally recorded on the 
Zhongdian Plateau of the Jinsha River drainage (upper Yangtze) 


Googl 


Figure 1 Bita Lake and the only native fish P. chungtiensis in it 





in China including Bita Lake, Napa Lake, Shudu Lake, the Naya, 
Xiaozhongdian and Geza Rivers! (Chen, 2010; Huang & Chen, 
1986). However, according to field surveys in recent years, it 
has been almost extirpated from its original localities, and only 
Bita Lake contains a relative large population (Jiang et al., 
2013). The fact that the P. chungtienensis is endemic to Bita 
Lake and is the only native species present makes this lake 
ecosystem unique. 


B 


A: the location of Bita Lake (from the screenshot of Google Earth); B: the aerial view of Bita Lake from the Northern hillside; C: the nearby view of Bita 


Lake from the Northern lakeshore; D: the living photo of P. chungtiensis. 


However, this special relationship of P. chungtienensis and 
Bita Lake has been interrupted by the introduction of non-native 
oriental weatherfishes species? at some point after 2009. The 
exotic species are generally threatening the native ones, 
particularly to those that have reduced populations with low 
genetic diversity (Magurran, 2009), and those in simple 





1 The population in Geza River was considered as a subspecies 
named P. chungtienensis gezaensis before with distinct morphological 
variations (Huang & Chen, 1986), but we now tentatively treated it as 
a distinct sister species P. gezaensis according to the regular practice 
in fish taxonomy that no more use of the category of subspecies. 

? These introduced oriental weatherfishes actually include two species 
with close relationships, similar morphology and ecological niches, 
Misgurnus anguillicaudatus and Paramisgurnus dabryanus. The two 
species were usually mixed in our samples so we simply refer to them 
as weatherfishes in this paper. 
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ecosystems with vulnerable trophic cascades (Estes et al., 
2011). The situation in Bita Lake is possibly one example of this. 
From 2008 to 2013, we carried out periodic surveys of the 
fishes in Bita Lake. This paper aims to provide a dynamic 
assessment of JP  chungtienensis and the introduced 
weatherfishes by successive monitoring data, and investigate 
the potential risks to this threatened species through its intrinsic 
genetic diversity and the extrinsic competition for food with 
weatherfishes. 


MATERIALS AND METHODS 


Acoustic surveys using scientific Echo sounders have 
increasingly become valuable and frequently-used methods in 
fish assemblages monitoring and fishery management 
(Fernandes, 2009). A portable Echo sounder HONDEX HE-51C 
operating at a frequency of 200 kHz was used during the 
surveys. Echo signals were recorded in the afternoon of each 


investigating day along the same route. Four non-overlapping 
transects (BT1, 2, 3, 4) were chosen to cover the different 
depths, crossing the relative shallow areas in the east (<3 m), 
the deep areas in the mid-south (>7 m), and most areas with an 
average depth of 4 m. Echo recordings were cumulated at each 
transect, and the averages of each transect at different times 
each year were calculated to produce the curve of annual 
fluctuations. The recordings of the four transects from different 
times were statistically compared using the one-way ANOVA 
test (SPSS package). Fish samples were captured at night at 
the same day in three sites of the lake body using traps, and 
three surrounding inlet creeks using a backpack electro-shocker. 
The numbers of P. chungtienensis and weatherfishes in each 
sampling were recorded to see the changes of fish composition. 
Most of the fishes were released the next morning, and only a 
few samples were fixed by formalin to do the length-weight 
relationship, age determination, sex ratio and gonad analyses 
(Jiang et al., 2013), and conduct genetic diversity and dietary 
composition analyses in this study. 

The genetic diversity of the mitochondrial cytochrome b (Cyt 
b) gene sequence of P. chungtienensis was analyzed. Total 
DNA of 81 specimens from Bita Lake were purified from 
alcohol-preserved fin by the standard methods (Sambrook et al., 
1989). The following cycling conditions were used in 
polymerase chain reaction (PCR): an initial denaturing step at 
95 °C for 4 min; 38 cycles of denaturing at 94 °C for 1 min, 
annealing at 52 °C for 1 min, and an extension at 72 °C for 1 min; 
with a final extension step of 72 °C for 10 min. Amplification 
primers of Cyt b were adopted from other Cyprinidae species 
as L14724 (5'-GACTTGAAAAACCACCGTTG-3’) and H15915 
(5'-CTCCGATCTCCGGATTACAAGAC-3’) (Xiao et al., 2001). 
All sequences have been deposited in GenBank (Accession No: 
KJ841795-KJ841875). Three sequences from GenBank 
(Accession No: FJ601043, AY463507, and AY463508, the last 
one from Xiaozhongdian River and the others from Bita Lake) 
were also downloaded and included in all of the analyses. One 
sequence of P. gezaensis (historically considered a subspecies) 
from the Geza River (GenBank Accession No. AY463506) was 
used as the outgroup to calculate the genetic distance and 
haplotype variations. Sequences were aligned with ClustalW 
(Thompson et al., 1994) and then calculated the pairwise 
distance (p-distance) in MEGA 6 (Tamura et al., 2013). The 
number of distinct haplotypes, haplotype diversity (h) (Nei, 1987) 
and nucleotide diversity (rr) (Nei & Tajima, 1981) were 
calculated using DnaSP version 5.0 (Librado & Rozas, 2009) to 
estimate the standard indices of genetic diversity. A Median- 
Joining (MJ) approach (Bandelt et al., 1999) was performed to 
visualize relationships among haplotypes in Network v.4.6.1.2 
(http:/Awww.fluxus-engineering.com/). 

Dietary composition of P.  chungtienensis and the 
weatherfishes was analyzed to compare their food niche 
relationships and degree of dietary overlap. Stomachs of all 
fishes collected were excised and fixed in a 10% formalin 
solution to avoid further digestion, and the total length (TL) and 
body weight (BW) were recorded. Prey items were identified to 
the lowest possible taxonomic level and enumerated based on 
their main physical characteristics under a microscope. All 


individuals and segments of each prey category were sorted 
dry and weighed jointly to the nearest 0.001 g. The percentage 
by number (%N), percentage by weight (%W), percentage 
frequency of occurrence (%F.O) and index of relative 
importance (IRI) were recorded to quantify the diet composition, 
and the overlap formula of Pianka (1973) and a modified PS 
index of Shorygin (1952) were used to assess the diet overlaps. 
These equations used to quantify the diet composition and 
assess the dietary overlap were listed below. 

Percentage by number (N%)=(number (n) of individuals of a 
prey category/total number (n) of individuals among all prey 
categories)x 100 (1) 

Percentage by weight (W%)=(weight (g) of individuals of a 
prey category/total weight (g) of individuals among all prey 
categories)x 100 (2) 

Percentage frequency of occurrence (F.O96)2(number (n) of 
stomachs containing a prey taxon/total number (n) of stomachs 
containing prey)x100 (Bowen, 1996; Hyslop, 1980) (3) 

IRI =%F.O (%N + %W) (Pinkas et al., 1971) (4) 
where IRI is the index of relative importance, %N, %W and 
%F.O are these defined at (1, 2, 3). 


Ox -Zra [$ Pi? Px? (Pianka, 1973) (5) 


where Oj is the overlap index of Pianka which indicated the 
niche overlap of species j and k, Pj and P; are the numerical 
proportions of the food item i used by the species j and k 
respectively. 


PS-Y (ab) 


where a and b are the percentages by number of the prey 
common to the two predators, and the percentage similarity (PS) is 
a summation of the smaller of the values of a and b for each prey. 
The index ranges from 0 (no overlap) to 100 (complete overlap). 


(Ivlev, 1961) (6) 


RESUTLS 


Population dynamics of P 
weatherfishes 

Eight effective Echo sounder surveys were made between 2009 
and 2012. The fluctuations of mean numbers each year in each 
transects can be seen in Figure 2A. Generally, the first two 
transects (BT1 & 2) recorded more individuals than that in the 
last two transects (BT3 & 4). It might be a result of the different 
characters of specific length and depth of each transect. These 
recordings of transects were not significantly different among 
different years (from 2009 to 2012, P=0.270), which shows no 
significant change in the quantity of fish in these years. A total of 
12 samples were used to compare the population dynamics of 
P. chungtienensis and weatherfishes can be seen in Figure 2B. 
The catch ratios of weatherfishes to P. chungtienensis varied 
from 0 (no weatherfishes) before August 2009 to a maximum of 
16.2 in August 2012. 


chungtienensis and 


Genetic diversity of P. chungtienensis 
A 1 140 base pair fragment of Cyt b sequence was amplified, 
with six variable sites and three parsimony-informative sites. 
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The overall p-distance within P. chungtienensis was 0.001, and 
the mean distance between P. chungtienensis and P. gezaensis 
was 0.017. Seven unique haplotypes were identified in all 84 
sequences. The h and rr values inferred from individuals of P. 
chungtienensis were 0.666 and 0.000 81, respectively. The 
Median-Joining network revealed a simple haplotype 
relationship for P. chungtienensis (Figure 3). There were at least 
20 mutations isolating P. chungtienensis and its sister species P. 
gezaensis. Haplotype 1 (H1) was a relatively ancestral 
haplotype of P. chungtienensis, because it was shared by 
individuals from Bita Lake and the Xiaozhongdian River. 
Another comparable haplotype, H2, was shared by 37 
individuals, and all other haplotypes radiated from H2 by only 
one or two mutations. 
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Figure 2 Field surveys of P. chungtiensis and weatherfishes in Bita 
Lake 

A: annual fluctuations in fish numbers (n) from Echo sounder recordings 
in four transects from 2009 to 2012; B: catch numbers (n) of P. 
chungtienensis and weatherfishes from 2008 to 2013. 


Diet composition and overlapping of P. chungtienensis 
and weatherfishes 

A total of 105 individuals of P. chungtienensis (90-221 mm in TL 
and 11.9-134 g in BW) were used to investigate stomach 
contents, of which 92 (87.696) individuals contained food. The 
sample of weatherfishes included 125 individuals (76-203 mm 
in TL and 3.8 to 74.5 g in BW) which 119 (95.296) had food. The 
diet composition (listed by N96, W%, F.O% and IRI) of P. 
chungtienensis and weatherfishes in Bita Lake is given in Table 
1. According to the IRI values for P. chungtienensis, algae 
(Spirogyra and Zygnema) were most important, followed in 
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order by Trichoptera (Philopotamus), Amphipoda (Gammarus 
bitaensis), and aquatic plants. These four items together 
represented 97.696 of the food intake. For weatherfishes, 
the sequence was Amphipoda (Gammarus bitaensis), algae 
(Spirogyra and Zygnema), and Diptera (larvae of Chironomus). 
These three items comprised 95.796 of the relative dietary 
importance. The overlap index of Pianka between P 
chungtienensis and weatherfishes was 0.69. According to a 
study of trophic relationships of fishes in an Ontario temperate 
lake (Keast, 1978), the O; value of 0.3 or less was insignificant, 
and one of 0.7 or more was considered high overlap of prey 
composition. The PS index of Shorygin between P 
chungtienensis and weatherfishes was 79.1696, which is much 
greater than the critical point of 60% (obvious similarity of prey 
composition (Blaber & Bulman, 1987)). Both the overlap index 
of Pianka and a modified PS index of Shorygin suggested a 
high degree of dietary overlap between P. chungtienensis and 
weatherfishes. 
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Figure 3 Median-joining network of P. chungtienensis based on 
seven haplotypes (H1 to H7) of Cyt b gene from 84 individuals in 
Bita Lake and the Xiaozhongdian River, with comparing to its 
sister species P. gezaensis in the Geza River 


DISCUSSION 


Non-native species have been introduced into new ecosystems 
primarily through human activity, either deliberately or 
unintentionally (Gozlan et al., 2010). The approach of 
deliberately introductions is mainly from societal demands for 
fish products for food aquaculture (51%), ornamental fish (21%), 
sport fishing (12%) and fisheries (7%) (Gozlan, 2008). Nearly 
8% of other introductions are accredited to accidental 
introductions, such as the escape from aquaculture installations, 
dispersal through ballast water, or introduced as a contaminant 
of a consignment of aquaculture species (Gozlan et al., 2010). 
The exotic weatherfishes presence in Bita Lake occurred due 
to an unusual pathway of biological introduction - the practice 
as “prayer animal release”. Bita Lake is a protected area in the 
Bitahai Nature Reserve, where intentional introduction of exotic 
species is strictly forbidden. However, the majority of citizens in 


Table 1 Diet composition of P. chungtienensis and weatherfishes in Bita Lake 





P. chungtienensis Weatherfishes 
Prey items N% W% F.O96 IRI N% W% F.0% IRI 
Amphipoda 33.9 51.02 50 4246 42.89 54.11 94.74 9 189.78 
Diptera 9.04 2.66 16.67 195.04 37.09 10.06 84.21 3970.5 
Trichoptera 53.67 11.78 83.3 5 451.99 1.22 1.01 5.26 11.73 
Copepoda — — — — 9.65 0.7 26.32 272.41 
Clcdocerans = = = — 4.57 0.35 31.58 155.37 
Gastropods — — — — 4.57 12.43 26.32 447.44 
Oligochaeta 1.67 1.93 33.3 119.88 — — — — 
Odonata 1.67 3.93 16.67 93.35 — — — — 
Algae 35 20.64 100 5564 45 21.34 100 6634 
Aquatic plants 16 8.04 66.67 1602.75 — — — — 


the area practice Tibetan Buddhism, when releasing animals 
into the wild are thought to generate good karma in 
metempsychosis. Buddhist monks have played an important 
role in promoting conservation efforts around the world (Chong, 
2012). However, religious belief can also lead to unsustainable 
practices without sound ecological knowledge (Gong et al., 
2012). The introduction of exotic weatherfishes into Bita Lake 
may be a negative example. Moreover, the active release of 
weatherfishes into the lake is sporadic and private, according to 
a conservation officer, as the believers usually use coat pockets 
or backpacks to transport the weatherfishes, and release them 
into the lake or the inflowing streams. 

Although the annual fluctuation from Echo recordings shows 
little or no significant difference in fish amount between 2009 
and 2012 in Bita Lake, the ratio of weatherfishes to P. 
chungtienensis by catch has been increasing since the former 
were first collected in August, 2009.  Ptychobarbus 
chungtienensis is an endemic fish with a current distribution 
confined to Bita Lake, and its genetic diversity revealed in this 
study is very low when even comparing the h and rr values to 
the pooled population of the flagship species of conservation, 
Chinese sturgeon Acipenser sinensis, one of the four Category 
| State protected fishes in China (h, 0.666 vs. 0.949; rr, 0.00081 
vs. 0.011) (Zhang et al., 2003). By contrast, weatherfishes 
thrive from the cold temperate to the subtropical regions in 
eastern Asia, and has been shown to expand its range very 
rapidly after introduction into new habitats throughout the world 
(Franch et al., 2008; United States Fish and Wildlife Service, 
2012). Another two lakes on the Zhongdian plateau, Shudu 
Lake and Napa Lake, were once populated by P 
chungtienensis. However, when various non-native fishes 
including weatherfishes were introduced into these lakes, the P. 
chungtienensis no longer occurs in them now. Now 
weatherfishes have been introduced into Bita Lake. Given their 
very successful survival strategy (United States Fish and 
Wildlife Service, 2012), and the high degree of dietary overlap 
with P. chungtienensis that we have demonstrated above, we 
are deeply concerned about the future existence of P 
chungtienensis. 


In order to further preserve P. chungtienensis from its last 
habitat, first of all, more effective measures are imperative to 
prevent other introduced species from Bita Lake, which need 
wisdom to balance the existing conflict between conservation 
and religious faith. Furthermore, continuous periodic monitoring 
and comparative studies are indispensable to find if there is any 
other critical factors that will put the P. chungtienensis in more 
danger. Last but not the least, some pre-arranged plans are 
also necessary before it's too late, such as the timely launching 
of artificial breeding research of P. chungtienensis. 
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ABSTRACT 


Fish of the superfamily Cobitoidea sensu stricto 
(namely loaches) exhibit extremely high diversity of 
color patterns, but so far little is known about their 
evolutionary mechanism. Melanocortin 1 receptor 
gene (MC1R) plays an important role during the 
synthesis of melanin and formation of animal body 
color patterns. In this study, we amplified and 
sequenced the partial MC1R gene for 44 loach 
individuals representing 31 species of four families. 
Phylogenetic analyses yielded a topology congruent 
with previous studies using multiple nuclear loci, 
showing that each of the four families was 
monophyletic with sister relationships of Botiidae+ 
(Cobitidae+(Balitoridae+Nemacheilidae)). Gene 
evolutionary analyses indicated that MC1R in 
loaches was under purifying selection pressure, with 
various sites having different dy/ds values. Both 
Botiidae and Cobitidae had lower dy/ds values than 
those of background lineages, suggesting their 
evolution might be strongly affected by purifying 
selection pressure. For Balitoridae and 
Nemacheilidae, both had larger dyds values than 
those of background lineages, suggesting they had a 
faster evolutionary rate under more relaxed selection 
pressure. Consequently we inferred that the 
relatively stable color patterns in Botiidae and 
Cobitidae might result from the strong purifying 
selection pressure on the MC1R gene, whereas the 
complicated and diverse color patterns in Balitoridae 
and Nemacheilidae might be associated with the 
relaxed selection pressure. Given the easy 
experimental procedure for the partial MC1R gene 
and its excellent performance in reconstructing 
phylogeny, we suggest this gene could be used as a 
good molecular marker for the phylogenetic study of 
fish species. 


Keywords: Gene evolution; Molecular phylogeny; 
Cobitoidea; Melanocortin 1 receptor gene 
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INTRODUCTION 


The superfamily Cobitoidea is a group of small- to medium- 
Sized benthic fish, composed of approximately 2896 of species 
of the order Cypriniformes, which is the largest group of 
freshwater fish in the world (Nelson et al., 2016). Depending on 
different authors, Cobitoidea includes variable families. Bohlen 
& Slechtová (2009) and Chen et al. (2009) congruently 
recognized the genus Ellopostoma as a distinct new family 
Ellopostomatidae, and proposed that Cobitoidea is composed 
of eight families (Catostomidae, Gyrinocheilidae, Botiidae, 
Vaillantellidae, Cobitidae, Ellopostomatidae, Nemacheilidae and 
Balitoridae). Kottelat (2012) raised genera Serpenticobitis and 
Barbucca to family rank, and established Serpenticobitidae and 
Barbuccidae. Both Serpenticobitis and Barbucca have been 
formerly included in Balitoridae (Bohlen & Slechtová, 2009; 
Slechtová et al., 2007), and consist of three and two species, 
respectively. While the most recent version of Fishes of the 
World (Nelson et al., 2016) follows Kottelat (2012), without 
more convincing evidence, we herein adopt the classification of 
Cobitoidea suggested by Slechtová et al. (2007), Bohlen & 
Slechtová (2009), and Chen et al. (2009). 

Cobitoidea sensu stricto, which excludes Gyrinocheilidae 
(algae eaters) and Catostomidae (suckers), is commonly known 
as the loaches (Nelson et al., 2016). This group of fish exhibit 
extremely high color pattern diversity, resulting in the observed 
abundant biodiversity of loaches. This makes them ideal 
subjects for studying adaptive evolution to environments, 
camouflage, kin recognition, and mate choice. Their beautiful 
color patterns also attract numerous fish enthusiasts, making 
them one of the most important groups of aquarium fish. A 
number of community websites about loaches have been 
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established by fish hobbyists, with the most detailed site 
(Loaches Online, http: //www.loaches.com/) containing 
comprehensive knowledge on their taxonomy, form, disease, 
and care. Additionally, some loach species are considered as 
important economic fish, including Misgurnus anguillicaudatus 
(mud loach), which is widely farmed in China (Zhang et al., 
2012). Therefore, loaches play significant roles in both 
evolutionary biology and social economy. To date, however, the 
phylogenetic relationships of loaches still remain controversial 
(Liu et al., 2012; Mayden et al., 2009; Slechtova et al., 2007; 
Tang et al., 2006). As for their color patterns, almost no 
investigations have been performed, and little is known about 
their evolutionary mechanisms. 

Among the many genes related to color patterns, 
melanocortin 1 receptor gene (MC1R) is one of the most widely 
studied, belonging to a family of G-protein-coupled receptors 
involved in various physiological processes in vertebrates and 
playing an important role in the synthesis of melanin and 
formation of animal body color patterns (Selz et al., 2007). The 
MC1R gene has been extensively investigated in numerous 
animal species, including reptiles (Corso et al., 2012; Cox et al., 
2013; Herczeg et al., 2010; Nunes et al., 2011; Rosenblum et 
al, 2004), birds (Baiao & Parker, 2012; Mundy, 2005), 
mammals (Ayoub et al., 2009; Lu & Zhang, 2001; Pérez et al., 
2013), and fish (Bar et al., 2013; Gross et al., 2009; Henning et 
al., 2010; Selz et al., 2007). Almost all fish species studied thus 
far have shown that the MC1R gene is conserved as a single 
copy and single exon protein-coding gene with a 966 bp long 
open reading frame (ORF), including seven transmembrane 
domains (Bar et al., 2013; Braasch et al., 2008; Henning et al., 
2010; Hofreiter & Schóneberg, 2010; Nunes et al., 2011), 
suggesting that this gene should be a good molecular marker 


Table 1 Samples used in the present study 


for reconstructing fish phylogeny. So far, however, only a few 
studies have applied this gene to mammal phylogenies (Ayoub 
et al., 2009; Pérez et al., 2013), and no report exists for fish 
phylogenies. The color patterns of loaches are highly diverse, 
which make them good subjects to study the evolution of the 
MC1R gene. In the present study, we amplified and sequenced 
the partial sequence of the MC7R gene in loaches, and 
analyzed their characteristics. Using this gene as a molecular 
marker, the phylogeny of loaches was reconstructed, and the 
evolution of the MC1R gene was then analyzed. 


MATERIALS AND METHODS 


Sample collection 

Because species of the families Ellopostomatidae and 
Vaillantellidae are distributed in Borneo, Sumatra, and Thailand, 
and only consist of two and three species, respectively, 
samples can be difficult to obtain. In this study, four families of 
loaches were involved. In total, 52 individuals representing 34 
species in Cypriniformes were selected for analyses of MC1R 
gene evolution, which included 44 individuals belonging to 31 
loach species and several outgroups, including one 
Gyrinocheilus aymonieri, one Hypophthalmichthys molitrix, and 
six Danio rerio individuals. To guarantee the accuracy of the 
targeted sequences, the cDNA sequences of four D. rerio 
individuals were downloaded from GenBank. In addition, we 
sequenced two D. rerio individuals, and all 46 other sequences 
in the present study. All samples were collected from the main 
drainages in China and aquarium markets, preserved in 95% 
ethanol, and deposited in the Institute of Hydrobiology, Chinese 
Academy of Sciences. Detailed information of all samples is 
listed in Table 1, with GenBank accession numbers included. 





Family Species and individuals Locality GenBank accession No. 
Botiidae Sinibotia reevesae Chishui, Guizhou Province KX120160 
Yasuhikotakia modesta Market in Wuhan KX120161 
Leptobotia elongata Market in Wuhan KX120164 
Leptobotia pellegrini Quzhou, Zhejiang Province KX120167 
Leptobotia microphthalma 1 Tucheng, Guizhou Province KX120162 
Leptobotia microphthalma 2 Tucheng, Guizhou Province KX120163 
Leptobotia taeniops 1 Yuanjiang, Hunan Province KX120165 
Leptobotia taeniops 2 Yuanjiang, Hunan Province KX120166 
Parabotia maculosus 1? Jianou, Fujian Province KX120168 
Parabotia maculosus 2° Jianou, Fujian Province KX120169 
Parabotia fasciatus Tucheng, Guizhou Province KX120170 
Cobitidae Cobitis sinensis 1 Wuyishan, Fujian Province KX120157 
Cobitis sinensis 2 Wuyishan, Fujian Province KX120158 
Pangio kuhlii Market in Wuhan KX120159 
Nemacheilidae Homatula wujiangensis Chishui, Guizhou Province KX120149 
Homatula variegata 4n Chishui, Guizhou Province KX120146 
Homatula variegata 2 Baoji, Shaanxi Province KX120148 
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Family Species and individuals 


Continued 





Nemacheilidae Homatula variegata 3° 
Schistura fasciolata 
Schistura incerta 
Traccatichthys taeniatus 1 
Traccatichthys taeniatus 2 
Triplophysa bleekeri 
Triplophysa sp. 1 
Triplophysa sp. 2 
Balitoridae Erromyzon kalotaenia 1° 
Erromyzon kalotaenia 2° 
Vanmanenia stenosoma 
Vanmanenia pingchowensis 
Vanmanenia lineata 
Vanmanenia hainanensis 
Vanmanenia caldwelli 1° 
Vanmanenia caldwelli 2° 
Formosania stigmata 
Formosania chenyiyui 1 
Formosania chenyiyui 2 
Pseudogastromyzon fasciatus 1 
Pseudogastromyzon fasciatus 2 
Pseudogastromyzon fangi 1° 
Pseudogastromyzon fangi 2° 
Beaufortia kweichowensis 
Jinshaia abbreviata ' 

Jinshaia sinensis ' 
Sinogastromyzon szechuanensis 
Outgroups Gyrinocheilus aymonieri 
Hypophthalmichthys molitrix 
Danio rerio 1 

Danio rerio 2 

Danio rerio 3 

Danio rerio 4 

Danio rerio 5 


Danio rerio 6 


Individuals with the same superscript letters share the same haplotype. 


DNA extraction, PCR amplification, and sequencing 

Total DNA was extracted from muscles following the salt- 
extraction procedure of Aljanabi & Martinez (1997), with some 
modification as per Tang et al. (2008). Five pairs of primers 
were newly designed according to the conserved regions of the 
aligned teleost MC1R gene sequences, and the lengths of the 


Locality GenBank Accession No. 
Baoji, Shaanxi Province KX120147 
Xichang, Sichuan Prov. KX120150 
Changting, Fujian Province KX120151 
Zhongshan, Guangzhou Province KX120156 
Zhongshan, Guangzhou Province KX120155 
Chishui, Guizhou Province KX120154 
Baoji, Shaanxi Province KX120152 
Baoji, Shaanxi Province KX120153 
Laibin, Guangxi Province KX120133 
Laibin, Guangxi Province KX120134 
Quzhou, Zhejiang Province KX120132 
Yuanjiang, Hunan Province KX120127 
Laibin, Guangxi Province KX120129 
Yuedong, Hainan Province KX120128 
Shaowu, Fujian Province KX120138 
Shaowu, Fujian Province KX120139 
Huangken, Fujian Province KX120142 
Changting, Fujian Province KX120130 
Changting, Fujian Province KX120131 
Guangze, Fujian Province KX120135 
Longyan, Fujian Province KX120136 
Laibin, Guangxi Province KX120140 
Longyan, Fujian Province KX120141 
Market in Wuhan KX120137 
Panzhihua, Sichuan Province KX120144 
Panzhihua, Sichuan Province KX120143 
Chishui, Guizhou Province KX120145 
Market in Wuhan KX120173 
Market in Wuhan KX120174 
Market in Wuhan KX120172 
Market in Wuhan KX120171 
AY161847* 
NM180970* 
BC162848* 
BC162836* 


*: sequences downloaded from GenBank. 


target segments ranged from 650 bp to 850 bp. However, only 
by using the primer set for the shortest target segment could we 
successfully obtain sequences for all samples. Thus, the primer 
set used was inMC1RF (5-AGCGTCAGYAAYGTGGTGGAGA- 
3) and inMC1RH (5-CGGTTCTGTACCTGCACAT-3?. 

The polymerase chain reaction (PCR) mixtures contained 3 uL of 
10x taq Buffer, 0.9 uL of dNTPs (10 mmol/L), 0.75 uL of each 
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primer (10 mmol/L), 0.15 uL of Taq polymerase, 1 uL of DNA 
template, and sterile distilled H2O to a final volume of 30 pL. 
The PCR thermal cycle profile was as follows: initial denaturation 
step at 94 °C for 5 min, followed by 35 cycles at 94 °C for 30 sec, 
56 °C for 45 sec, 72 °C for 1 min, and a final extension at 72 °C 
for 8 min. The amplified fragments were purified and sequenced 
by Sangon Biotech (Shanghai) Co., Ltd. To confirm the 
accuracy of sequencing, sequencing reactions were performed 
from both ends of each fragment using the same primers as 
PCR amplification. 


Sequence variation and phylogenetic analyses 

Multiple alignments of sequences were performed using 
ClustalX 2.1 (Larkin et al., 2007; Thompson et al., 1997) 
alignment editor with the default settings, with all aligned 
sequences then translated into amino acid residues in MEGA 
6.0 (Tamura et al., 2013) to test for sequencing mistakes. The 
alignments were verified by eye in SEAVIEW (Galtier et al., 
1996). Saturation of genes and codon positions were assessed 
using the ‘transitions and transversions vs. divergence’ graphic 
function in DAMBE 5 (Xia, 2013). Nucleotide compositions and 
pairwise genetic distances were calculated based on the 
Kimura 2-parameter model with standard errors estimated 
using 1 000 bootstrap replicates. Molecular phylogenetic 
relationships were estimated using Bayesian inference (BI) and 
maximum likelihood (ML) methods. The optimal model of 
nucleotide evolution for BI analyses was identified using the 
Bayesian Information Criterion (BIC) as estimated in jModelTest 
2.1.7 (Darriba et al., 2012; Guindon & Gascuel, 2002). The 
TPM1uf+G model was selected as the best model, so we set 
the parameter nst=6 with rates=gamma when running the BI 
analyses. 

Bayesian inference (Bl) was carried out using MrBayes 3.1.2 
(Ronquist & Huelsenbeck, 2003). Two independent analyses 
with four simultaneous Markov chains were run for 4 000 000 
generations, sampling every 1 000 generations, with a total of 4 
001 trees each. We discarded all samples obtained during the 
first 1 000 000 generations as "burn-in", and a 5096 majority- 
rule consensus tree with posterior probability values for each 
node was obtained from the remaining 3 001 trees (Figure 1). 

Maximum likelihood (ML) analysis was conducted using 
RAxML 8.1.13 (Stamatakis, 2014). We searched for the best 
scoring ML tree using a general time-reversible nucleotide 
model (GTR) with gamma-distributed rate variation among sites 
(G) and invariable sites (Il). A rapid bootstrap analysis with 1 
000 replicates was used to assess the relative robustness of 
node support. 


Analyses of MC1R gene evolution 

The w ratio (defined as the ratio of non-synonymous vs. 
synonymous substitution rates, i.e., w=dy/ds) is a measure of 
natural selection acting on a protein. Simplistically, values of 
w<1, 21, and >1 means negative purifying selection, neutral 
evolution, and positive selection, respectively (Yang, 2007). Site 
models and branch models were separately adopted to 
estimate whether the MC1R gene in the present study had 
experienced different selection pressures along different sites 


284  www.zoores.ac.cn 


and different lineages. The CODEML program within the PAML 
package was used to assess parameters in models of 
sequence evolution and to test relevant hypotheses (Yang, 
2007). In the following analyses, the outgroups were 
excluded, and the input tree was an unrooted tree with 
trifurcation at the root. 

For site models, we examined three pairwise codon-based 
substitution models to assess w values for all codon sites: MO 
(one-ratio) vs. M3 (discrete w), M1a (nearly neutral) vs. M2a 
(positive selection), and M7 (8 distribution) vs. M8 (B distribution 
and a fraction of sites with w>1). Likelihood ratio tests (LRTs) 
were performed to compare the fit of two pairwise models. It is 
assumed that twice the log likelihood difference between 
nested models (2AInL) follows a X^ distribution with the number 
of degrees of freedom equal to the difference in the number of 
free parameters (Whelan & Goldman, 1999). When LRTs 
indicate positive selection, the Bayesian empirical Bayesian 
(BEB) method (Yang et al., 2005) can be used to calculate the 
posterior probabilities of positively selected codons. 

For branch models, based on the phylogenetic tree yielded in 
the present study (Figure 1), we divided loaches into four 
lineages corresponding to four families (Balitoridae, 
Nemacheilidae, Cobitidae, and Botiidae). The following steps 
were performed: (1) w values of the MC1R gene dataset were 
estimated for all lineages using a one-ratio model (all lineages 
share a common w value); (2) each lineage was separately 
treated as a foreground lineage compared with the remaining 
lineages (background lineages), and  two-ratio models 
(foreground lineage and background lineages have different w 
values) were independently conducted for each lineage; (3) four 
pairs of LRTs were examined between the one-ratio and two- 
ratio models for four lineages. 


RESULTS 


Characteristics of MC1R gene sequences 

After alignment, 591 bp of the MC1R gene sequences was 
used for analysis. No indels or deletions were detected. In total, 
43 haplotypes were identified from 52 individuals. For the 44 
loach individuals, 38 haplotypes were recognized. Except for 
Jinshaia abbreviata and J. sinensis sharing the same haplotype, 
all other shared haplotypes occurred among different 
individuals of the same species. The average nucleotide 
composition for all species was A=19.5%, T=28.5%, G=22.1%, 
and C=29.9%. The content of A+T (48.0%) was lower than that 
of C+G (52.0%). Strong compositional biases against A existed 
at the third position (only 7.6%). For the 591 bp sequences, 226 
sites were variable, of which 194 were parsimony informative. 
The average ratio of transitions/transversions (Ti/Tv) was 2.084. 
Plots of the absolute transitions and transversions versus TrN 
distance showed that the transitions and transversions for all 
three codon positions did not reach saturation (not shown). 


Phylogenetic relationships of Cobitoidea 

Maximum likelihood and BI analyses yielded almost the same 
topologies, except that the lineage (representing one botiid 
subfamily, Botiinae) formed by Sinibotia reevesae and 
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Figure 1 ML tree estimated using RaxML for all individuals based on partial MC1R gene sequences 
Values at the nodes correspond to support values for ML/BI methods (only values >50 or 0.50 shown). *: a support value of 100 or 1.00; -: values of «50 
or 0.50. Species names are followed by individual codes as in Table 1. 
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Yasuhikotakia modesta had a different position. In the ML tree, 
this lineage clustered with botiid subfamily Leptobotiinae, 
formed by Leptobotia and Parabotia species, and together they 
formed the family Botiidae. Whereas in the BI tree, Botiinae was 
located in the basal position of all loaches, but with a very low 
posterior probability (only 0.44). The ML tree with node support 
values of both methods labeled is shown in Figure 1. 

The phylogenetic tree indicated that each of the four loach 
families was monophyletic and had high support values 
(bootstrap value (BP) higher than 94 in ML, and all posterior 
probability (PP)=1.00 in Bl), except for Botiidae. The sister 
relationships among the four families was Botiidae* (Cobitidaet+ 
(Balitoridae+Nemacheilidae)). Balitoridae and Nemacheilidae 


formed a sister group with medium support values (BP=91 in ML, 
and PP=0.89 in BI), located at the crown of the phylogenetic 
tree. The family Cobitidae was at an intermediate position. 
Within Botiidae, Botiinae was distantly related to the 
Leptobotiinae with a genetic distance of 0.092+0.013. For 
balitorid species, subfamilies Balitorinae and Gastromyzontinae 
were both monophyletic and formed monophyletic Balitoridae. 
Within Gastromyzontinae, however, the species of genera 
Vanmanenia, Formosania, and Pseudogastromyzon were not 
monophyletic. The genetic distance between these two 
subfamilies was not very high (0.073+40.010). The pairwise 
mean genetic distances of the MC1R gene among the loach 
groups are listed in Table 2. 


Table 2 Pairwise mean genetic distance (lower triangle) and standard error (upper triangle) based on the Kimura 2-parameter model for 


loach groups and outgroups 








Group code 
Family or subfamily Group code 

1 2 3 4 5 6 7 
Gastromyzontinae 1 0.010 0.012 0.012 0.014 0.013 0.018 
Balitorinae 2 0.073 0.011 0.012 0.014 0.012 0.017 
Nemacheilidae 3 0.114 0.093 0.013 0.015 0.012 0.017 
Cobitidae 4 0.114 0.090 0.121 0.013 0.012 0.016 
Botiinae 5 0.123 0.120 0.142 0.118 0.011 0.016 
Leptobotiinae 6 0.116 0.093 0.104 0.101 0.092 0.015 
Outgroups 7 0.190 0.174 0.186 0.175 0.167 0.150 


Bold indicates the genetic distance between two subfamilies of the same family. 


Evolution of the MC1R gene in Cobitoidea 

Using the topology yielded by the ML method, we analyzed the 
evolution of MC1R in Cobitoidea. Analysis based on the one- 
ratio model MO showed that the average w value of the MC1R 
gene in all loaches was 0.034 9, much less than 1, suggesting 
that the average selection pressure on this gene was negative 
or purifying. Using LRT analyses, comparisons of the pairwise 
models indicated that the M3 model (having discrete w values 
among sites; -InL=3 295.024 728) was significantly better than 
the MO model (having the same w value of all sites; -/nL=3 
366.540 405) (P=0.000) at estimating w values among sites; 


however, there was no obvious difference between M1a 
(nearly neutral; -/nL-3 301.622 991) and M2a (positive 
selection; -/nL-3 301.622 991) models (P=1.000), suggesting 
that no amino acid site was obviously under positive 
selection, although different sites had variable 
evolutionary rates. However, comparison also indicated 
that the M8 model (-/nL-3 297.026 100) was significantly 
better than the M7 model (-/nL-3 304.204 399) (P=0.001), 
and the BEB procedure identified amino acid sites 12Q, 20Q, 
and 176H as possibly being under positive selection, though 
each posterior probability (Pr) was less than 95% (59.2%, 
73.096, and 87.396, respectively) (Table 3). 


Table 3 Likelihood ratio tests (LRTs) between selected CODEML codon substitution models for the MC1R gene of Cobitoidea 


2AInL 
143.033 2 


Compared models InL 
M3-MO M3: -3 295.024 728 


MO: -3 366.540 405 


M2a-M1a M2a: -3 301.622 991 0 
Mia: -3 301.622 991 
M8-M7 M8: -3 297.026 100 14.356 6 


MT: -3 304.204 399 


For variation of the w values among the loach lineages, 
pairwise comparisons of LRTs between so called "foreground" 
and "background" lineages showed that two pairs were 
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df P value Positive selected sites (posterior probability) 
4 0.000 

2 1.000 

2 0.001 12Q (59.2%), 20Q (73.0%), 


176H (87.3%) 


detected with significant variations. The mean w value of 
Balitoridae (0.055 1) was significantly higher than that of the 
other three families (denoted as non-balitorids, 0.026 0) 


(P=0.006), whereas the mean w value of Cobitidae (0.002 8) 
was significantly lower than that of the non-cobitids (0.040 5) 
(P=0.00). For Botiidae, its mean w value was lower (0.025 2) 
than that of the non-botiids (0.039 8) (P=0.096). The mean w 


value of the nemacheilids (0.046 0) was slightly higher than 
that of the non-nemacheilids (0.032 1) (P=0.211). The 
pairwise comparison results for the four lineages are listed in 
Table 4. 


Table 4 LRTs performed to detect heterogeneous selection regimes among lineages for the MC1R gene 








Branch-specific model Parameter estimates df InL 2AInL P value 

All lineages 

MO (one-ratio) w=0.034 9 87  -3 366.540 405 

Balitoridae* 

two-ratio vs. one-ratio Wo=0.026 0 w1=0.055 1 88 | -3 362.719 906 7.640 998 0.006<0.01** 


Nemacheilidae? 
two-ratio vs. one-ratio Wo=0.032 1 w1=0.046 0 
Cobitidae? 

two-ratio vs. one-ratio Wo=0.040 5 w1=0.002 8 
Botiidae? 


two-ratio vs. one-ratio (970.039 8 (470.025 2 


88 | -3365.759 399 1.562 012 0.211 


88 . —-3358.153 163 16.774 48 0.00«0.01** 


88 | -3365.156 736 2.767 338 0.096 


a: foreground lineage; w4: value of dw/ds for foreground lineages; wo: value of dw/ds for background lineages. 


DISCUSSION 


Evolutionary variation of MC1R among sites and loach 
lineages 

Presently, dozens of studies have indicated that the MC1R 
gene is highly conserved among vertebrates, and variations are 
closely associated with pigmentation differences (Hubbard et al., 
2010; Mundy, 2005; Nunes et al., 2011; Rosenblum et al., 2004), 
with only a few exceptions (Cheviron et al., 2006; Dorn et al., 
2011). Such findings are important in clarifying the origin of new 
species or local adaptation within species (Hubbard et al., 
2010). Compared with other vertebrates, teleost fish have more 
diverse color patterns (Parichy, 2003) Among numerous 
cypriniform fish, loaches are one of the most important groups 
with colorful pigment patterns. 

The present results showed that in loaches the MC1R gene 
was under negative selection pressure overall (average value 
of w=0.034 9«1). Analyses based on site models indicated that 
the M3 model was significantly better than the MO model 
(P=0.000), suggesting that various sites of the MC1R gene had 
disparate selection pressure. However, no amino acid site was 
obviously under positive selection since no significant difference 
existed between the Mia (nearly neutral) and M2a (positive 
selection) models (P=1.000). In contrast to the M1a-M2a 
comparison, the results of the M7-M8 comparison were 
significant and model M8 detected three positive amino acid 
sites, though with low posterior probabilities. Yang (2007) 
determined that the w estimate under M8 might often produce 
false positives and the M1a-M2a comparison could be more 
robust. Therefore, given the results of the M1a-M2a comparison, 
evolution of the MC1R gene in loaches does not appear to be 
affected by positive selection. 

As for different loach lineages, analyses based on branch 
models showed that both Cobitidae and Botiidae had much 


lower w values than the background lineages, with Cobitidae 
reaching a level of significant difference and Botiidae showing 
no significant difference. Since differences in w values reflect 
differences in the level of constraint (Rosenblum et al., 2004), 
the lower w values suggest that MC1R gene evolution in both 
Cobitidae and Botiidae might be strongly affected by purifying 
selection pressure. For Balitoridae and Nemacheilidae, 
however, both have higher w values than those of the 
background lineages, with the former at a significant level. This 
suggests that both families might have experienced a more 
relaxed selection pressure and relatively faster evolutionary rate 
than that of the background lineages. Among the four loach 
families, the botiid species exhibit the most regular color pattern, 
commonly with black vertical bars on the body sides. Species of 
Cobitidae have several types of color patterns, though typically 
show Gambetta pigment lines, as displayed in the genus 
Cobitis, which are defined as five color zones from the dorsal to 
lateral sides (Chen & Chen, 2005; Chen et al., 2015). For 
balitorid and nemacheilid species, the color patterns are 
extremely complicated and diverse, and include vertical bars, 
horizontal bands, spherical blotches, scattered spots, and 
irregular cloud-like patches. Consequently, we inferred that the 
relatively stable color pattern in Botiidae and Cobitidae might 
partially result from strong purifying selection pressure on the 
MC1R gene, whereas the diverse color patterns in Balitoridae 
and Nemacheilidae might be associated with relaxed selection 
pressure on this gene. 


Phylogenetic application of the MC1R gene 

As mentioned above, in the fish species analyzed so far, the 
MC1R gene is reported to be a conserved single copy and 
single exon protein-coding gene with a total length of more than 
800 bp, making it a good molecular marker for phylogenetic 
research (Li et al, 2007). Until this study, however, no 
phylogenetic application of this gene has been reported for fish. 
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In our research, partial MC1R gene sequences (591 bp) were 
adopted for establishing the phylogeny of Cobitoidea, which 
yielded a well-supported phylogenetic tree, except for Botiidae 
with low support value. The present topology was consistent 
with previous studies using multiple nuclear loci (Chen et al., 
2009; Liu et al., 2012; Mayden et al., 2009), and showed that 
Botiidae included two distantly related clades, corresponding to 
the two subfamilies Leptobotiinae and Botiinae proposed by 
Slechtova et al. (2006). The genetic distance between these 
two subfamilies was 0.092+0.013, higher than that between the 
two subfamilies of Balitoridae (0.073+0.010). Therefore, it is 
reasonable to divide Botiidae into these two subfamilies. Within 
Balitoridae, the three genera, Vanmanenia, Formosania and 
Pseudogastromyzon, were not monophyletic. In our previous 
study (Liu et al., 2012), the monophyly of Vanmanenia and 
Formosania could not be supported either. Further study is 
needed to confirm the phylogeny of these three genera using 
more species samples, molecular markers, and morphological 
characters. 

It is worth noting that the partial segment of the MC1R gene 
used in this study was easy to amplify using genomic DNA as a 
template with general PCR procedures, and easy to sequence. 
Therefore, this gene segment could be used as a good 
molecular marker for the phylogenetic study of fish species in 
the future. 
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ABSTRACT 


Triplophysa daochengensis sp. nov. is described 
from the Daocheng River, a northern tributary of the 
Jinsha River in Sichuan Province, China. The new 
species can be distinguished from its congeners by 
the following characters: body smooth and scales 
absent; lateral line complete; caudal peduncle 
compressed, depth unchanging; head length equal 
to caudal-peduncle length; lower jaw shovel-shaped; 
dorsal-fin origin anterior to pelvic-fin origin and 
closer to the tip of the snout than to the caudal-fin 
base, last unbranched ray hard; pelvic-fin tip not 
reaching anus; posterior chamber of gas bladder 
absent; intestine of spiral type with three winding 
coils. 


Keywords: Nemacheilidae; Triplophysa; New species; 
Sichuan 


INTRODUCTION 


The genus Triplophysa is wide spread on the Qinghai-Tibet 
plateau and in adjacent areas. In total, 127 species of 
Triplophysa have been reported, with 109 occurring in 
China (Yan et al., 2015). Triplophysa can be distinguished 
from other Nemacheilidae genera by the character 
combination of close together nostrils, and sexual 
dimorphism in which males have patches of tubercles on 
both sides of the head and an agglomeration of tubercles 
on the dorsal surfaces of the branched and broadened 
pectoral-fin rays (Yang et al., 2011). 

The Jinsha River from Yushu in Qinghai to Yibin in Sichuan 
marks the upper reaches of the Yangtze River. Twenty species 
of Triplophysa have been recorded from the Jinsha River basin so 
far: Triplophysa aquaecaeruleae, T. anterodorsalis, T. bleekeri, T. 
brevibarba, T. brevicauda, T. daqiaoensis, T. grahami, T. 
hialmari, T. leptosoma, T. markehenensis, T. ninglangensis, T. 
orientalis, T. pseudostenura, T. stenura, T. stolickai, T. venusta, 
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T. xichangensis, T. xigiensis, T. yajiangensis and T. yaopeizhii 
(Ding, 1993; Ding & Lai, 1996; Guo et al., 2012; He et al., 2012; 
Prokofiev, 2001; Wu & Wu, 1992; Xu & Zhang, 1996; Yan et al., 
2015; Zhu, 1989). 

From March 1995 to June 2015, we collected 2 959 
specimens of nemacheilid loaches in the Jinsha River basin. 
Comprehensive study of these specimens and detailed 
comparison with species previously recorded from the Jinsha 
River (Ding, 1993; Ding & Lai, 1996; Guo et al., 2012; He et al., 
2012; Prokofiev, 2001; Wu & Wu, 1992; Xu & Zhang, 1996; Yan 
et al., 2015; Zhu, 1989) unveiled a new species, which is 
described herein. 


MATERIALS AND METHODS 


Specimens were captured using hand nets and fish traps. The 
geographic coordinates of the specimen collection sites were 
obtained using a GPS Garmin eTrex handheld device. 
Specimens were fixed in a 1096 formalin solution in the field, 
and after 5 d were transferred into a 596 formalin solution for 
long-term storage. Measurement methods followed Wu & Wu 
(1992) and Prokofiev (2007) and were made with digital calipers 
recorded to its nearest 0.1 mm. Lateral head length was from 
snout tip to the most posterior point of the opercle. 

Abbreviations used in this paper are: SL, standard length; HL, 
lateral head length; CWNU, College of Life Sciences, China 
West Normal University, Nanchong, China; KIZ, Kunming 
Institute of Zoology, Chinese Academy of Sciences, Kunming, 
China; IHB, Institute of Hydrobiology, Chinese Academy of 
Sciences, Wuhan, China. 
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RESULTS 
Triplophysa daochengensis sp. nov. (Figures 1) 


Holotype: CWNU 201406001, male, 85.0 mm SL; Daocheng 
River, at Shengmu Village, Daocheng County, Sichuan 
Province, China; N29°02'55", E101°30'45"; 3 540 m a.s.l.; 
collected by Yan-Shu Guo, Jun Yang, and Ming-Hao Luo on 
June 16, 2014. 


Paratypes: 21 specimens, 69.0-98.0 mm SL. CWNU 


E 


Figure 1 Triplophysa daochengensis sp. nov. 


201006001-006 collected by Yan-Shu Guo on June 01, 2010; 
CWNU 201406002, KIZ 2014005963, IHB 201406004, CWNU 
201406005-201406006, KIZ 2014005964, CWNU 201406008- 
201406009 collected by Yan-Shu Guo, Jun Yang, and Ming- 
Hao Luo on June 16, 2014; same collection site as the holotype. 
IHB 201406010, CWNU 201406011-201406014, KIZ 
2014005965, IHB 201406016 collected by Yan-Shu Guo, Jun 
Yang, and Ming-Hao Luo on June 16, 2014; the collection site 
was a wetland by the Daocheng River, at Chahua Village, 
Daocheng County, Sichuan Province, China; N29°02'04", 
E100°23'43"; 3 700 m a.s.l. 





A, B, D, E, F, G: holotype, CWNU 201406001, 85.0 mm SL, male (A: lateral view, life coloration; B: ventral view; D: dorsal view; E: ventral view of head; 
F: lateral view of head; G: pectoral-fin); C: paratypes, CWNU 201406012, 92.0 mm SL, female, lateral view; H: paratype, CWNU 201006006, 96.8 mm 
SL, male, air bladder; | and J: paratype, CWNU 201406008, 95.5 mm SL, female (I: ventral view of intestines; J: dorsal view of intestines). 


Diagnosis: Body smooth and scales absent; lateral line 
complete; caudal peduncle compressed, depth unchanging; head 
length equal to caudal-peduncle length; lower jaw shovel-shaped; 
dorsal-fin origin anterior to pelvic-fin origin and closer to the tip of 
the snout than to the caudal-fin base, last unbranched ray hard; 
pelvic-fin tip not reaching anus; posterior chamber of gas bladder 
absent; intestine of spiral type with three winding coils. 

Description: Morphometric and meristic characteristics are 
provided in Table 1. Body elongate, caudal peduncle 
compressed and depth to caudal fin direction unchanging. Body 
surface smooth and scaleless. Lateral line complete with 86-95 
pores. Vertebrae 4+40 (n=2). 

Snout blunt, with length equal to postorbital head length. 
Anterior and posterior nostrils located adjacently and close to 
anterior rim of orbital. Mouth inferior. Lips thick and furrowed; 
upper lip without a median incision; lower lip with v-shaped 
central notch; lower lip surface with shallow furrows. Lower jaw 
shovel-shaped, with a sharp edge, not covered by lower lip. 
Inner rostral barbel not extending to corner of mouth, outer 
rostral barbel extending to a vertical through nostril, maxillary 
barbel reaching a vertical through posterior margin of eye. 

Dorsal fin with a concave distal margin, last unbranched ray 
hard. Dorsal-fin origin anterior to pelvic-fin origin and closer to 
the tip of the snout than to the caudal-fin base. Pectoral-fin tip 


exceeds the midpoint between pectoral and pelvic-fin origins. 
Pelvic-fin tip not reaching anus. Anal fin with straight distal 
margin. Caudal-fin deeply emarginate, lower lobe slightly longer 
than the upper one. 

Gas bladder with an anterior chamber fully enclosed in bony 
capsule; posterior chamber absent. Stomach "u'"-shaped. 
Intestine winding spiral type with three coils. 

Color pattern in life: Color of the dorsum and flanks yellowish- 
brown, abdomen white. Eight to twelve black brown transverse 
saddles on back, 10-16 black spots on the side of the body 
post-medianly. All fins grayish yellow. Two to three transverse 
rows of small black spots on the dorsal fin and 2-4 rows of small 
dark spots on the caudal fin. 

Color pattern of preserved specimens: Dorsum and flanks 
gray, abdomen whitish beige. Eight to twelve black transverse 
saddles on dorsum and 10-16 black spots on the side of the 
body post-medianly. 

Sexual dimorphism: Upper and lower patch of tubercles in 
pre- and suborbital areas in males; upper patch is a short strip 
skewed upwards located sub-antero-orbitally; lower patch is 
triangle-shaped and located at the base of the maxillary barbel 
in front of the operculum. In males, the first to sixth branched 
pectoral-fin rays possess a patch of tubercles dorsally; pectoral- 
fin is shorter but wider than in females. 
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Table 1 Main morphometric and meristic characters of T. daochengensis sp. nov., T. brevibarba, T. stolickai and T. yajiangensis 


Triplophysa daochengensis sp. nov. 





Paratypes n=21 


T. brevibarba n-8 


T. stolickai n=6 


T. yajiangensis n-22 





Holotype 





Characters Range MeantSD Range MeanzSD Range MeanzSD Range MeanzSD 
Total length (mm) 101.0 81.0- 80.0-102.0 68.0-78.5 80.0-125.0 
117.0 
SL (mm) 85.0 69.0-98.0 67.0-87.0 59.0-65.5 69.0-107.0 
Percentage of SL (96) 
Head length 20.5 19.3-22.5 21.040.9 20.4-23.0 21.541.0  22.9-23.7 23.440.5  20.8-228 21.9+0.7 
Body depth 12.4 11.0-14.5 12.1£1.3  20.1-22.4 = 21.2:1.0  11.9-16.0 14.142.1  13.0-16.5 14.3112 
Body width 10.6 10.4-14.5 11.341.6 . 17.1-20.3 19.0£1.4  11.0-16.0 13.5425  112-16.0 12.841.5 
Pectoral-fin length 16.5 16.0-19.8 18.241.2 . 18.8-20.9 19.940.9 18.6-19.1 18.8403 16.8-19.1 17.8+0.8 
Pelvic-fin length 15.9 14.2-16.7 15.840.8  14.4-16.4 15.440.9 14.4-16.0 15.2408 144-169 15.8+0.9 
Caudal-peduncle length 20.5 19.3-22.5 20.840.9 16.1-20.9 18.842.7  23.7-26.3 24.641.5  24.3-25. 24.7+0.5 
Predorsal length 48.2 48.1-49.4 48.8402  51.5-55.3 54441.9 50.8-53.4 51.741.5 43.5-48.3  46.2+1.5 
Distance pectoral to pelvic-fin 29.4 27.2-29.9 28.6:0.7 33.8-36.6  35.941.2 31.44-35.41 33.241.9  25.6-29.3 27.041.3 
Percentage of HL (%) 
Head depth 51.4 50.0-53.9 51.941.1 58.8-66.7 63.5425  48.1-50.0 48.8+1.0 50.0-56.3 52.542.4 
Head width 62.8 57.5-64.7 61.6+2.3 77.1-80.6 78.141.8 51.9-58.1 55.7:3.3  60.0-72.7 63.6+5.0 
Snout length 40.0 39.0-42.5  41.240.7 46.7-51.4 48.742.0 35.7-40.5 38.2:3.3  39.9-41.8 41.140.7 
Postorbital head length 40.0 39.0-42.5 41.240.7 36.1-41.2 39.641.9 39.5-46.4 43.2440  40.2-41.9 41.2+0.6 
Inner rostral-barbel length 17.1 15.0-18.4  17.141.6 16.7-20.0 17.641.2  17.9-19.4 18.640.8 21.7-26.7 23.942.4 
Outer rostral-barbel length 28.6 20.0-29.4 24.843.1 22.2-26.7 = 23.541.7  17.9-22.2 19.842.2  30.0-36.44  33.0£2.1 
Maxillary-barbel length 28.6 20.0-29.4 27.043.6  16.7-26.7 17.641.2  21.4-29.0 25.443.8 31.0-39.5 34.6+3.6 
Pectoral-fin length/Distance 56.0 58.6-68.4 63.8448  52.4-62.0 56.2435 54.0-59.5 56.6+2.8 58.1-69.2 66.8+5.0 
pectoral to pelvic-fin (%) 
Head length/Caudal-peduncle 100.0 94.4- 99.140.7 = 107.1-128.6 113.2412.2 87.1-100.0 95.74+7.4 81.5-92.3 88.943.7 
length (%) 100.0 
Caudal-peduncle depth/length (%) 40.0 36.1-42.1 39.942.9 42.9-46.4  42.843.0 22.6-25.8 24.541.7 32.5-40.8  36.6£3.0 
Dorsal-fin rays ii-7 iii-7 iii-8 iii-7 iii-8 
Anal-fin rays ii-5 ii-5 ii-5 ii-5 ii-5 
Pectoral-fin rays i -10 i -10 i-11-12 i-9 i-10 
Pelvic-fin rays i-7 i-7 i-7 i-7 i-8 
Gill rakers 15 15-17 21-23 16-18 13-15 
Distribution: Currently known in Daocheng River, a northern 
tributary of the Jinsha River in Sichuan Province, and its DISCUSSION 


associated wetlands (Figure 2). 

Ecology: Inhabits areas close to river banks and wetlands 
with slow water current bottom composed of gravel and mud; 
submerged waterweeds present (Figure 3). Algae and aquatic 
invertebrates were found in the stomachs of six dissected 
specimens. Several other species were collected at the same 
time: T. stolickai, T. brevicauda, Misgurnus anguillicaudatus, 
Schizothorax kozlovi and Schizopygopsis malacanthus. 

Etymology: The specific name is derived from its collection 
location in Daocheng County. 
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Triplophysa daochengensis can be distinguished from T. 
orientalis, T. xichangensis, T. venusía, T. obscura, T. 
pseudoscleroptera, and T. scleroptera by lower jaw shovel- 
shaped vs. spoon-shaped and posterior chamber of the gas 
bladder absent vs. present (Ding, 1994; Wu & Wu, 1992; Zhu, 
1989; Zhu & Wu, 1981). 

Triplophysa daochengensis can be distinguished from T. 
aquaecaeruleae, T. bleekeri, T. brevicauda, T. dagiaoensis, T. 
grahami, T. leptosoma, T. ninglangensis, T. pseudostenura, 


Iaa eyssure 


Figure 2 Collection site of Triplophysa daochengensis sp. nov. 





e Specimen collection site 


Scale: 1 : 1200 000 





Figure 3 Habitat of Triplophysa daochengensis sp. nov. 


Mainstem and backwaters of the Daocheng River, Daocheng County, photograph taken on June 16, 2014. 


T. xigiensis, T. yaopeizhii, T. angeli, T. pappenheimi, T. robusta, 
and T. siluroides by lower jaw shovel-shaped vs. spoon-shaped 
and intestines of winding spiral type with three coils vs. one 
zigzag loop (Ding, 1993, 1994; Ding & Lai, 1996; Guo et al., 
2012; He et al., 2012; Prokofiev, 2001; Wu & Wu, 1992; Xu & 
Zhang, 1996; Zhu, 1989). 

Triplophysa daochengensis, T. alexandrae, T. hialmari, T. 
stenura and T. lixianensis share a shovel-shaped lower jaw, but 
the caudal-peduncles of 7. alexandrae, T. hialmari, T. stenura, 
and T. lixianensis are round and taper obviously vs. 
compressed and depth unchanging (Figure 1, Figures 4A, H) 
(Chen et al., 2004; Ding, 1994; He et al., 2008; Prokofiev, 2001; 
Wu & Wu, 1992; Zhu, 1989). 

Triplophysa daochengensis and T. anterodorsalis share a 
shovel-shaped lower jaw, but can be distinguished by last 


unbranched dorsal-fin ray hard vs. soft; caudal-fin deeply 
emarginate vs. emarginated; intestine of spiral type with three 
winding coils vs. one zigzag loop; pectoral-fin length/distance 
from pectoral to pelvic-fin 58.6%-68.4% vs. 68.7%-89.8%; and 
postorbital head length equal to snout length vs. shorter than 
snout length (Figure 1, Figure 4D) (Ding, 1994; Yan et al., 2015; 
Zhu, 1989). 

Triplophysa stolickai, T. crassilabris and T. daochengensis 
share a shovel-shaped lower jaw; however, T. stolickai and T. 
crassilabris can be distinguished from T. daochengensis by the 
last unbranched dorsal-fin ray soft vs. hard; dorsal-fin origin 
further to the tip of the snout than to the caudal-fin base vs. vice 
versa; pelvic-fin tip reaching or exceeding the anus vs. not 
reaching the anus (Table 1, Figure 1, Figures 4E, G) (Chen et 
al., 2004; Ding, 1994; Wu & Wu, 1992; Zhu, 1989). 
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Figure 4 Eight related species of Triplophysa daochengensis sp. nov. 


A: T. stenura (CWNU 201407006, 98.5 mm SL, male); B: T. brevibarba (CWNU 201112079, 92.5 mm SL, male); C: T. markehenensis (CWNU 
201108021, 83.2 mm SL, male); D: T. anterodorsalis (CWNU 201112066, 69.2 mm SL, male); E: T. stolickai (CWNU 201406069, 65.5 mm SL, male); F: 
T. yajiangensis (CWNU 201407001, 107.0 mm SL, female). G: T. crassilabris (CWNU 200807036, 68.5 mm SL, male); H: T. lixianensis (CWNU 


201108032, 103.5 mm SL, male). 


Triplophysa daochengensis and T. markehenensis share a 
shovel-shaped lower jaw, but can be distinguished by the 
posterior chamber of the gas bladder being absent vs. present; 
dorsal-fin origin nearer to the tip of the snout than to the caudal- 
fin base vs. vice versa; caudal-peduncle depth to caudal fin 
direction unchanged vs. slightly decreasing; and pelvic-fin tip 
not reaching the anus vs. reaching or exceeding the anus 
(Figure 1, Figure 4C) (Wu & Wu, 1992; Yan et al., 2015; Zhu, 
1989; Zhu & Wu, 1981). 

Triplophysa daochengensis can be distinguished from T. 
brevibarba by dorsal-fin origin nearer to the tip of the snout than 
to the caudal-fin base vs. vice versa; rakers on the inner side of 
the first branchial arch 15-17 vs. 21-23; and caudal-peduncle 
depth towards caudal-fin unchanged vs. obviously tapering 
(Table 1, Figure 1, Figure 4D) (Ding, 1993; Yan et al., 2015). 

Triplophysa daochengensis can be distinguished from T. 
yajiangensis by caudal-peduncle depth to caudal-fin direction 
unchanged vs. slightly decreasing; pelvic-fin tip not reaching the 
anus vs. reaching or exceeding the anus; and head length 
equal to caudal-peduncle length vs. shorter than caudal- 


peduncle length (Table 1, Figure 1, Figure 4F) (Yan et al., 2015). 


Triplophysa daochengensis can be distinguished from T. 
nujiangensa by snout length equal to postorbital head length vs. 
shorter than postorbital head length; predorsal length 48.1%- 
49.4% of SL vs. 71.1%-77.2%; last unbranched dorsal-fin ray 
hard vs. soft; pectoral-fin length/distance from pectoral to pelvic- 
fin 58.6%-68.4% vs. 47.5%-63.5%; and caudal-peduncle depth 
to caudal-fin direction unchanged vs. slightly decreased (Chen 
et al., 2004). 

Triplophysa daochengensis can be distinguished from T. 
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tanggulaensis by head length equal to caudal-peduncle length 
vs. longer than caudal peduncle; last unbranched dorsal-fin ray 
hard vs. soft; pelvic-fin tip not reaching the anus vs. exceeding 
the anus; vertebrae 4+40 vs. 4+37-38; and rakers on the inner 
side of the first branchial arch 15-17 vs. 10-13 (Zhu, 1989; 
Chen et al., 2004). 
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CWNU 201108032-052, CWNU 201407067-083, Zagunao 
River, Lixian County, Aba Prefecture. T. markehenensis. 
(topotype) CWNU 201108001-026, Maerke River of the upper 
Dudu River, Banma County, Qinghai Province. T. obscura 
CWNU 198706001-024, CWNU 199508001-016, Bilonjiang 


(Cypriniformes: Cobitidae). Acta Zootaxonomica Sinica, 21(3): 377-379. (in 
Chinese) 


Yan SL, Sun ZY, Guo YS. 2015. A new species of Triplophysa Rendahl 
(Cypriniformes, Nemacheilidae) from Sichuan Province, China. Zoological 
Research, 36(5): 299-304. 


Zhu SQ, Wu YF. 1981. A new species and a new subspecies of loaches of 
the genus Nemachilus from Qinghai Province. Acta Zootaxonomica Sinica, 
6(2): 221-224. (in Chinese) 
Zhu SQ. 1989. The Loaches of the Subfamily Nemacheilinae in China 
(Cypriniformes, Cobitidae). Nanjing: Jiangsu Science and Technology 
Press, 69-129. (in Chinese) 


River, Jiangzha Xiang, Roergai County; CWNU 199208001-036, 
Yellow River basin, Xiaman Xiang, Roergai County. T. 
orientalis CWNU 201108047-055, CWNU 201407061-075, 
Yalong River basin, Shiqu County. T. pappenheimi CWNU 
200308001-012, Yellow River, Tangke Xiang, Roergai County; 
CWNU 200608001-008, Yellow River basin, Hoyuan County, 
Aba Prefecture. T. pseudoscleroptera CWNU 200208013-031, 
Yellow River basin, Tangke Xiang, Roergai County; CWNU 
200807001-023, Yellow River basin, Axi Xiang, Roergai County. 
T. pseudostenura (topotype) CWNU 199507051-052, CWNU 
201007013-029, Yalong River, Ganzi County, Ganzi Prefecture. 
T. robusta CWNU 198706025-069, Bilonjiang River, Dolie 
Xiang, Roergai County. T. scleroptera CWNU 200608001-019, 
Yellow River basin, Hoyuan County. T. siluroides CWNU 
201008017-020, Yellow River, Aba County, Aba Prefecture; 
Roergai Yellow River 201008032-056, Yellow River, Tangke 
Xiang, Roergai County. T. stenura CWNU 201407001-016, 
Jinsha River basin, Shiqu County, Ganzi Prefecture; CWNU 
201309001-006, Jinsha River basin, Yuanpu Xiang, Jiangda 
County, Tibet Autonomous Region, China. T. stolickai CWNU 
201407001-004, Jinsha River basin, Shiqu County; CWNU 
201406068-077, Daocheng River, Daocheng County T. 
xichangensis (topotype) CWNU 200812048-56, CWNU 
201509001-026, Anning River, Mianning County. T. xiqiensis 
(topotype) CWNU 201112052-076, Xiqi River, Jinsha River 
basin, Zhaojue County, Liangshan Prefecture. T. yajiangensis 
CWNU 201009005 (holotype), CWNU 20100901-03 (paratype), 
Yalong River, Yajiang County. T. yaopeizhii (topotype) CWNU 
201309007-019, Jinsha River basin, Yuanpu Xiang, Jiangda 
County, Tibet Autonomous Region; CWNU 201309001-005, 
Jinsha Jiang basin, Gonjue County, Tibet Autonomous Region. 
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ABSTRACT 


In September and October 2015, a new species of 
the family Nemacheilidae, Triplophysa tianxingensis 
sp. nov., was discovered from underground water in 
Qiubei County, Yunnan Province, China. It can be 
distinguished from all other troglobiotic Triplophysa 
species occurring in Yunnan by the following 
combination of characters: eyes small, a little 
degenerated; barbels longer; ventral profiles greatly 
convex; pectoral fin short, attaining a third of the 
distance from the pectoral-fin base to pelvic fin base; 
body with many brown blotches; caudal peduncle 
with fin fold; caudal fin shallowly forked, and free 
posterior chamber of swim bladder cylindrical. 


Keywords: Triplophysa; Cave fish; New species; 
Yunnan 


INTRODUCTION 


The genus Triplophysa Rendahl, 1933, one of the largest in the 
family Nemacheilidae, comprises 170 currently recognized 
species, with most occurring in the Qinghai-Tibet Plateau and 
its adjacent regions (Froese & Pauly, 2016; Zhu, 1989). The 
genus Triplophysa is distinct from other genera by the following 
combination of characters: close together nostrils, dumbbell- 
shaped bony bladder, bony posterolateral bladder capsule, and 
marked sexual dimorphism. In males, the lateral sixth to eighth 
branched rays of the pectoral fins are thickened, broadened 
and covered by breeding tubercles; breeding tubercles are also 
present on the sides of the head, extending from the eye almost 
to the insertion of the maxillary barbels (Yang, 1990; Zhu, 1989). 
Species divergence is very high in numerous isolated rivers and 
lakes. Due to their small body size and little interest to the 
fishery industry, their taxonomy has lagged compared with other 
larger-sized fish, and many undescribed species have been 
discovered in recent years. This genus occurs widely in central 
Asia, from Afghanistan and Baluchistan through the High Asian 
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Region to the Balkhash and Uvs-Nuur lakes, Outer Mongolia 
and China. In China, in addition to the Qinghai-Tibetan plateau 
and Inner Mongolia, the genus also occurs in Sichuan, Yunnan, 
Chongqing, Hunan, Guizhou and Guangxi provinces. To date, 
only seven troglobiotic Triplophysa species have been 
described in the Yunnan Province, namely, T. gejiuensis Chu & 
Chen (1979), T. nanpanjiangensis Zhu & Cao (1988), T. 
Shilinensis Chen et al. (1992), T. aluensis Li & Zhu (2000), T. 
yunnanensis Yang (1990), T. xiangshuigingensis Li (2004), and 
T. qiubeiensis Li et al. (2008). Among them, six species 
demonstrate degenerated small eyes or no eyes at all. 

In September and October 2015, we collected specimens of 
the genus Triplophysa, representing an undescribed species, 
from a seasonal pond connected to subterranean waters in 
Longtao Village (altitude: 1 257 m), Tianxing Township, Qiubei 
County, Yunnan Province, China (Figure 1). We herein provide 
a description of the new species and comparison to other 
troglobiotic Triplophysa species occurring in Yunnan. 


MATERIALS AND METHODS 


Specimens were collected with a cage net and were fixed in 
10% formalin and then transferred to 75% ethanol for long-term 
storage. Fish were cataloged and stored in the collections of the 
Heilongtan Reservoir of Shilin County (20151002001-6, 
20151004001, 20151004002) and the Kunming Institute of Zoology, 
Chinese Academy of Sciences (2015005390, 2015005391). 
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Figure 1 Sampling site of Triplophysa tianxingensis (marked with *) 


Measurements were taken point to point with digital calipers 
and data were recorded to 0.1 mm. Counts and measurements 
were made on the left side of specimens whenever possible. 
Methods of counts and measurements followed Kottelat (1990) 
and Chu & Chen (1989). Subunits of the head are percentages of 
head length (HL). Head length itself and measurements of body 
parts are expressed as proportions of standard length (SL). 


RESULTS 


Triplophysa tianxingensis sp. nov. (Figures 2, 3) 

Holotype 

20151002001, 118.0 mm SL, 136.0 mm total length; seasonal 
pond connected to subterranean waters in Longtao Village, 
N23°56'29.0"; E104°13'04.0"(altitude: 1 257 m), Tianxing 
Township, Qiubei County, Yunnan Province, China. 

Paratypes 

20151002002-6, 20151004001-2, (2015005390, 2015005391), 





Figure 2 Triplophysa tianxingensis sp. nov. (20151002001, 118.0 mm SL) 
A: Lateral view; B: Ventral view; C: Dorsal view. 
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Figure 3 Head of Triplophysa tianxingensis sp. nov. (20151002001, 
118.0 mm SL) 

A: Lateral view; B: Ventral view; C: Dorsal view. 

70.0-114.0 mm SL, 81.0-132.0 mm total length; nine specimens, 
collected with holotype. 
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Diagnosis 

This species can be distinguished from all other troglobiotic 
Triplophysa species occurring in Yunnan by the following 
combination of characters: eyes small, a little degenerated; 
barbels longer; ventral profiles greatly convex; pectoral fin short, 
attaining one third of the pelvic fin base; body with many brown 
blotches; caudal peduncle with fin fold; caudal fin shallowly 
forked, free posterior chamber of swim bladder cylindrical. 


Description 
Dorsal fin iii, 8; anal fin ii, 5; pectoral fin i, 9; pelvic fin i, 5; caudal fin 
with 16 branched rays; 2+12~13 gill rakers on the outside of the 
first gill arch (two specimens); 4+38 vertebrae (two specimens). 
Proportional measurements of Triplophysa tianxingensis sp. 
nov. are shown in Table 1. Body elongated; cylindrical anteriorly 
and a little compressed posteriorly; dorsal profile slightly convex 
and ventral surface distinctly convex. Dorsal-fin origin near 
midpoint of body, somewhat closer to caudal-fin base than to 


snout tip. Head slightly depressed, profile of head triangular, 
width greater than height at nape; head length greater than 
body depth. Snout obtuse and depressed, length less than 
postorbital length of head. Nose valve triangular. Both anterior 
and posterior nares close together, the latter larger than the 
anterior one and enclosing anterior posteriorly; anterior nostril 
pierced in the front side of a flattened tube; its tip reaching the 
posterior margin of the posterior nostril. Interorbital space 
convex. Mouth inferior and arched. Both lips thick and well- 
developed, with a few very shallow furrows. Lower lip broad, 
with a median notch. No processus dentiformis presenting on 
upper jaw, no corresponding notch on lower jaw. Eyes small 
and not visible in ventral view. Inner rostral barbel reaching 
vertical of the anterior margin of the nostril; outer rostral barbel 
surpassing vertical of the posterior margin of the eye. Maxillary 
barbel surpassing vertical of the posterior margin of the eye. 
Bony capsule of the airbladder dumbbell-shaped, posterior 
process cylindrical. Intestine simple, with no coil. Peritoneum 
light yellow ventrally and yellow gray dorsally. 


Table 1 Proportional measurements of Triplophysa tianxingensis sp. nov. 


2015100 2015100 2015100 2015100 2015100 2015100 2015100 2015100 2015005 2015005 





ID'numger 2001 2002 2003 2004 2005 2006 4001 4002 390 391 Range Mean 
Total length 136.0 120.0 1260 1320 1230 1300 1210 1090 810 1000 810-1360 117.8 
Standard length (SL) 118.0 105.0 1110 1140 1080 1120 1060 930 700 860 700-1180 1023 
Percentage of SL 
Head length (HL) (36) 2290 219 216 219 213 214 217 237 214 233 213233 224 
Body depth at dorsal fin 

= 169 143 162 167 185 161 170 183 157 174 143-185 16.7 
origin (%) 


Predorsal length (%) 55.9% 57.1 52.3 54.4 53.7 
Caudal peduncle length (96) 18.6 19.0 19.8 19.3 19.4 
Caudal peduncle depth (%) 9.3 10.5 9.9 10.5 10.2 


Percentage of HL (%) 


Snout length 44.4 50.0 43.8 46.0 43.5 
Eye diameter 5.6 6.5 6.3 6.0 4.3 
Interorbital width 22.2 30.4 20.8 24.0 21.7 


Percentage of caudal peduncle length 


Caudal peduncle depth (%) 50.0 55.0 50.0 54.5 52.4 


Distal margin of dorsal fin truncate. Pectoral fin extending 
about one third of the distance between the pectoral and pelvic 
fin origins. Axillary pelvic lobe present as vestige, indistinct. 
Pelvic fin reaching near half the distance between the pelvic 
and anal fin origins. Tip of the pelvic fin not attaining anus. Anus 
situated in front of anal fin origin, in the middle distance 
between the posterior tip of the pelvic fin and anal fin origin. 
Caudal fin shallowly forked, tips round. 

Skin smooth, no scales on the whole body. Lateral line 
complete, straight. 


Sexual dimorphism 
No sexual dimorphism was observed in the ten specimens. 
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52.7 53.8 57.0 52.9 57.0 52.3-57.0 54.7 


17.0 20.8 17.2 20.0 18.6 17.0-20.8 19.0 


8.9 9.4 9.7 10.0 10.5 8.9-10.5 9.9 


41.7 43.5 40.9 40.0 45.0 40.0-50.0 43.9 


4.2 4.3 4.5 6.7 5.0 4.2-6.7 5.3 


20.8 17.4 18.2 23.3 20.0 17.4-24.0 21.9 


52.6 45.5 56.3 50.0 56.3 50.0-56.3 52.3 


Usually, sexual dimorphism can be seen in mature individuals 
of Triplophysa species, namely, the pre-suborbital area might be 
thickened, and the pectoral fin might become broadened, 
thickened, stiffened and covered by blunt tubercles. These 
phenomena were not present in the ten individuals of T. 
tianxingensis sp. nov., which could be due to the sampling time 
being outside the propagation season of the species. 

Color pattern 

Color in life. Ground color of body light yellow, slightly lighter 
ventrally. Body with many brown blotches. Dorsal fin dark gray, 
with no distinct transverse gray bars. Caudal fin with three 
indistinct brownish transverse bars. Pectoral fin gray dorsally 
and pale ventrally. Pelvic and anal fin pale. Rostral barbel 


brown to yellow. Maxillary barbel light yellow (Figure 4). 





Figure 4 Live specimens of Triplophysa tianxingensis sp. nov. 


Color in formalin. Ground color of body light yellow, slightly 
lighter or pale ventrally. Body with many brown blotches. Dorsal 
fin dark grey, with no distinct transverse bars. Caudal fin with 
three indistinct dark brownish transverse bars. Pectoral fin gray. 
Pelvic and anal fin light grey. Rostral barbel dark grey. Maxillary 
grey (Figure 2). 


Habitat description 

This species was collected from a seasonal pond connected 
to subterranean waters in Longtao village, N23°56'29.0"; 
E104°13'04.0"(altitude: 1 257 m), Tianxing Township, Qiubei 
County, Yunnan Province, China. Water temperature 22 °C, 13: 
00 pm, Oct, 2, 2015. This region is famous for its well- 
developed karst landforms, as well as seasonal and 
subterranean rivers. These special geological conditions might 
provide very limited habitat for each species in this area (Figure 
5A, B). 


Distribution 

So far, this species is known only from the type locality 
(Figure 5A, B). This seasonal pond is connected with 
subterranean waters belonging to the Qingshuijiang River basin 
(a tributary of the Nanpangjiang River, upper reaches of the 
Pearl River). 


Etymology 

The specific name, tianxingensis (K & ), is based on the 
Tianxing Township sampling site, but is also an allusion to the 
unique color pattern of T. tianxingensis sp. nov., specifically, the 
striking brown blotches resembling stars (tianxing means stars 
in the sky in Chinese). 


DISCUSSION 


By 2008, seven troglobiotic Triplophysa species had been 
described in Yunnan Province, namely, T. gejiuensis Chu & 
Chen (1979), T. nanpanjiangensis Zhu & Cao (1988), T. 
shilinensis Chen et al. (1992), T. aluensis Li & Zhu (2000), T. 
yunnanensis Yang (1990), T. xiangshuigingensis Li (2004), and 





Figure 5 Habitat of Triplophysa tianxingensis sp. nov. 
A: Seasonal pool; B: Dry cave through which the seasonal pool 
connects with subterranean waters. 


T. qiubeiensis Li et al. (2008). Among them, six demonstrate 
degenerated small eyes or no eyes at all. The new species 
described in this paper can be differentiated from T. 
qiubeiensis, T. gejiuensis and T. shilinensis by the character 
of degenerated, but distinct eyes; differentiated from T. 
nanpanjiangensis, T. aluensis, T. yunnanensis and T. 
xiangshuiqingensis by the profile of head triangular and 
caudal peduncle with fin fold. 

An identification key is summarized below based on the of 
external characters of the eight troglobiotic Triplophysa species 
occurring in Yunnan. 


Key to eight troglobiotic Triplophysa species occurring in 
Yunnan: 
1(6) Eyes completely degenerated 
2(3) Head conical, nose valve triangular (Qiubei, Yunnan)..... 
C———Á— À T. qiubeiensis Li et Yang 
3(2) Head a little flat, nose valve beared 
4(5) Pelvic fin i, 5; anus situated immediately in front of anal 


fin origin (Gejiu, Yunnan) .................... T. gejiuensis Chu et Chen 
5(4) Pelvic fin i, 6; length between anus and anal fin origin 
about 1-2 mm (Shilin, Yunnan)......... T. shilinensis Chen et Yang 


6(1) Eyes normal or a little degenerated, small 
7(8) Eyes normal (Zhanyi, Yunnan)............................... sess 
——————— T. nanpanjiangensis Cao et Zhu 


Zoological Research 37(5): 296-300, 2016 299 


8(7) Eyes a little degenerated, small 
9(14) Head depressed or a little depressed, caudal peduncle 
with no fin fold 

10(11) Anterior part of body depressed (Luxi, Yunnan)........... 
EE T. aluensis Li et Zhu 

11(10) Anterior part of body cylindrical 

12(13) Posterior part of body with small scales (Yiliang, 
Yuman) asrar T. yunnanensis Yang 


13(12) No scales on body (Shilin, Yunnan).............................. 
EE ETE tre ERE En T. xiangshuiqingensis Li 
14(9) Profile of head triangular, caudal peduncle with fin fold 
(Qiubei, Yunnanin r a e 
A EETA T. tianxingensis Li, Yang et Chen, sp. nov. 


Comparisons of the characters among the eight troglobiotic 
Triplophysa species occurring in Yunnan are presented in Table 2. 


Table 2 Comparisons of characters among the eight troglobiotic Triplophysa species occurring in Yunnan 


T. nanpanjiangensis T. gejiuensis T. yunnanensis T. shilinensis T. aluensis 


T. xiangshuigingensis T. qiubeiensis T. tianxingensis 





Specimens 5 8 5 2 1 1 11 10 
Sampling site Zhanyi Gejiu Yiliang Shilin Luxi Shilin Qiubei Qiubei 
Dorsal fin iii, 7-8 iii, 7-8 iii, 7 iii, 7 iii, 7 iii, 6 iii, 7 iii, 8 
Anal fin ii, 5 iii, 4-6 iii, 5 iii, 5 iii, 5 iii, 5 iii, 5 ii, 5 
Pectoral fin i, 9-10 i, 10-11 i, 9-10 i, 8-10 i, 9 i, 9 i, 7-9 i, 9 
Ventral fin i, 6 i, 5 i, 7 i, 6 i, 6 i, 6 i,5 i, 5 
Caudal-fin rays 16 14-15 15-16 14 13 14 14-15 16 
Scales Absent Absent Present Absent Absent Absent Absent Absent 
Head shape Depressed Depressed Depressed Depressed Depressed Depressed Conical Triangular 
Anterior part of Cylindrical Cylindrical Cylindrical Cylindrical Depressed Cylindrical Cylindrical Cylindrical 
body 
Eyes Normal Completely Small Completely Degenerated Small Completely Small 
degenerated degenerated degenerated 
Caudal fin Emarginate Forked Emarginate Forked Forked Forked Forked Shallowly forked 
Distal margin Truncate Truncate Emarginate Truncate Truncate Emarginate Emarginate Truncate 
of dorsal fin 
Bar/Blotch Present Absent Present Absent Present Present Absent Present 
Position of Closer to Closer to At midpoint Closer to At midpoint Closer to At midpoint Closer to caudal 
dorsal fin origin snout tip caudal fin of body caudal fin of body snout tip of body fin base 
base base 
Fin fold Absent Absent Absent Absent Absent Absent Absent Present 
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ABSTRACT 


Piwi-interacting RNA (piRNA) plays an important role 
in the gonadal development and maintenance of 
Teleostei. In this study, piRNA libraries derived from 
the adult gonads of Japanese flounder (Paralichthys 
olivaceus) were generated using next-generation 
sequencing technology. Using zebrafish piRNAs as 
a reference, 5 865 unique candidate piRNAs were 
identified; 289 candidate piRNA clusters (PRCs) 
were generated from the above piRNAs. Among the 
isolated candidate PRCs, a total of 38 ovary-specific, 
45 ovary-bias, 24 testis-specific, and 131 testis-bias 
PRCs were found. The relative expression levels of 
seven PRCs were validated through quantitative 
reverse transcription-polymerase chain reaction. The 
results of this study will help facilitate exploration of 
the development and maintenance of the phenotypic 
sex mechanism in P. olivaceus. 


Keywords: Paralichthys olivaceus; Ovary; Testis; 
piRNA 


INTRODUCTION 


In mammals, the corresponding protein of piwi-interacting RNA 
(piRNA) interacts with a class of germ cell-restricted small 
RNAs (26-31 nucleotides in length) in the gonads. piRNAs 
show distinctive localization patterns in the genome and are 
predominantly grouped into 20—90-kilobase clusters, whereas 
the long stretches of small RNAs are derived from only one 
strand. In humans and rats, with major clusters occurring in 
syntenic locations, both the abundance of piRNAs in germline 
cells and male sterility of Piwi mutants suggest the possible role 
of piRNAs in gametogenesis (Aravin et al., 2006; Girard et al., 
2006). In zebrafish (Danio rerio), the piwi gene is expressed in 
both the testis and ovary and is also a component of a germline 
structure called "nuage". The loss of piwi function results in a 
progressive loss of germ cells due to apoptosis during larval 
development. Many of these small RNAs are derived from 
transposons; therefore, piRNAs are possibly involved in the 


Science Press 


silencing of repetitive elements in gonads (Houwing et al., 
2007). In general, genome-wide piRNA expression and function 
have been examined in fruit flies (Drosophilidae) (Aravin et al., 
2007), humans (Homo sapiens) (Sigurdsson et al., 2012), mice 
(Mus musculus) (Zheng et al., 2010), zebrafish (Houwing et al., 
2007; Huang et al, 2011), and planaria (Prosorhochmus 
claparedii) (Palakodeti et al., 2008). Furthermore, piRNAs are 
found in clusters throughout the genome, and each cluster may 
contain as few as 10 or up to thousands of piRNAs. The 
clustering of piRNAs is highly conserved across species, 
whereas the sequences are not (Malone et al., 2009). 

The gender of Japanese flounder (Paralichthys olivaceus) is 
easily altered by water temperature during the sex 
determination stage (Wen et al., 2014). Female P. olivaceus 
grow faster than males; thus, investigations on the mechanisms 
of sex differentiation can be potentially useful in increasing 
production of farmed P. olivaceus. Study on sex manipulation in 
P. olivaceus has been a topic of research for some time. Sex- 
related genes have been studied to elucidate the sexual 
molecular mechanism in flounder, and include the cyp19a gene, 
which was first cloned and found to be predominantly 
expressed in female flounders (Kitano et al., 1999); some 
transcription factors, such as foxl2 (Yamaguchi et al., 2007), 
dmrt1 (Wen et al., 2014), and dmrí4 (Wen et al., 2009), and 
other sex-related genes, such as cyp17 (Ma et al., 2012). 
However, due to the limited data on these genes, information 
regarding their expression and regulation mechanisms is also 
incomplete. As a representative of non-coding RNA, 
approximately 20.0% and 13.1% of isolated microRNAs 
(miRNAs) are preferentially expressed in the testis and ovary, 
respectively, which indicates that miRNAs might play a key role 
in the regulation of gene expression in P. olivaceus gonads (Gu 
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et al., 2014). However, the differential expression and function 
of piRNAs in fish sex development and maintenance remain 
poorly investigated. Identifying the expression patterns of 
piRNAs in P. olivaceus gonads is the first step in clarifying their 
function in gonadal development. Thus, this study aimed to 
identify candidate piRNAs in sexually mature gonads of P 
olivaceus through next-generation sequencing, which provides 
not only sequences of low-abundance piRNAs, but also 
quantitative data on the frequency of sequencing reads that 
correspond to the abundance of piRNAs in the ovary and testis. 


MATERIALS AND METHODS 


Fish maintenance 

One-year-old P. olivaceus fish with average individual body 
weights of 450—500 g were obtained from the Beidaihe Central 
Aquaculture Experiment Station, Chinese Academy of Fishery 
Sciences. All fish were acclimatized to indoor recirculation 
systems at the Shanghai Ocean University for two weeks. The 
gonads were examined through histological sectioning before 
RNA extraction, and only healthy fish were chosen for RNA 
isolation. All experiments were performed according to the 
Experimental Animal Management Law of China and as 
approved by the Animal Ethics Committee of Shanghai Ocean 
University. 


Histology and RNA isolation 

All fish were killed by overdose with tricaine methanesulfonate 
(MS-222, 200 mg/L). Mature ovaries and testes were collected 
separately from five females and five males. For histology, the 
gonads of females and males were fixed with Bouin's solution, 
then embedded in paraffin (Sigma Aldrich, http: 
Iwww.sigmaaldrich.com). Cross sections were cut at 5 uim and 
stained with hematoxylin and eosin. Tissues were isolated, 
immediately frozen in liquid nitrogen, and stored at —80?C. Total 
RNA was extracted using a mirVana'" miRNA isolation kit 
(Ambion; http:  //www.thermofisher.com/cn/zh/home/brands/ 
invitrogen/ambion.html) in accordance with the manufacturer’s 
instructions. Each testis or ovary was separately subjected to 
RNA extraction. Five RNA samples (1.5 ug per sample) derived 
from the same type of gonad were pooled for library 
construction and subsequent sequencing. The quantity and 
purity of total RNA were monitored using Nanodrop ND-1000 
(Nanodrop Technologies; http: //www.thermoscientific.com/content/ 
tfs/en/product/nanodrop-lite-spectrophotometer.html) and denaturing 
gel electrophoresis. The total RNA was stored at —80?C. 


Small RNA library construction and sequencing 

Two small complementary DNA (cDNA) libraries were 
generated from the mixture of total RNAs of mature testes 
and ovaries, respectively. A combination of RNA samples 
was used to avoid possible expression bias of piRNAs in 
different individuals. A small RNA fraction (15-31 nt) was 
isolated using 15% polyacrylamide/8 mol/L urea/0.5xTris/borate/ 
ethylenediaminetetraacetic acid gel electrophoresis and ligated 
with proprietary 5' (GTTCAGAGTTCTACAGTCCGACGATO) 
and 3' (TCGTATGCCGTCTTCTGCTTGT) adaptors using T4 
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RNA ligase. The short RNAs were then converted to cDNA 
through reverse transcription polymerase chain reaction (RT- 
PCR). The results were verified through Pico green staining and 
fluorospectrophotometry and then quantified using Agilent 2100; 
afterward, the cDNA libraries were sequenced using HiSeq2000 
(Illumina; http: //www.illumina.com/) in accordance with the 
manufacturer's recommended protocols for small RNA 
sequencing (NEXTflex"" Small RNA-seq kits; BIOO Scientific 
Corp; http: //www.biooscientific.com/). 


Bioinformatics analysis of sequencing data 

Raw sequences were filtered to remove low-quality reads 
(sequences with more than 1096 sequencing errors), and those 
that passed the quality filter were trimmed to remove the 
adaptor sequences. Sequences with 15-31 nt in the two 
combined libraries were then selected as raw small RNA 
sequences. The raw small RNA sequences were mapped to the 
piRNA database (Lakshmi & Agrawal, 2008) with zebrafish 
piRNAs used as the templates. The 5 865 candidate piRNAs 
are listed in Supplemental Table 1. Highly conserved unique 
piRNAs with different lengths were placed in the same piRNA 
clusters (PRCs), and total reads of all unique piRNAs were 
considered as the reads of the PRCs. The expression levels of 
the PRCs were estimated using RPKM (reads per kilobase of 
exon model per million mapped reads). Differential expression 
analyses of the two samples were performed using the DESeq 
R package (Anders et al., 2013). P-values were adjusted using 
the Benjamini and Hochberg method. A corrected P-value of 
0.05 was set as the threshold for significantly differential 
expressions. The P-values were also adjusted using Q-values; 
Q-value «0.01 and log; (fold change) >1 were set as the 
threshold for significantly differential expressions by default 
(Storey, 2003). 


Table 1 Primers used for piRNA qRT-PCR analysis 


Primer name Sequence(5’-3’) 





piR_0048450 GATGAACCGAACGCCGGGTTAA 
piR_0006070 AAGATGGTGAACTATGCCTGGG 
piR_0019030 GCAAATCGGTCGTCCGACCTGG 
piR_0015958 GATGAACCGAACGCCGGGTTAA 
piR_0055443 TATCATCCAGGTACACAAAGACA 
piR_0056629 TCCATTTTCAGGGCTAGTTGATT 
piR_0050867 TGATGATAGGTGTGAGTGCTACA 
Po5Ss CCATACCACCCTGAACAC 

Po5Sa CGGTCTCCCATCCAAGTA 


Quantitative RT-PCR analysis of piRNAs 

The expression profiles of nine PRCs were investigated with 
real time RT-PCR using a miScript reverse transcription kit and 
miScript SYBR Green PCR kit (Qiagen: http: //www. 
transmedchina.com/) in accordance with the manufacturers' 
instructions. Quantitative RT-PCR (qRT-PCR) was analyzed 
using iCycler (Bio-Rad; http: //www.bio-rad.com/). The 20 uL 
reaction mixture consisted of 10 uL of SYBR Premix Taq (2x), 


0.4 uL of piRNA-specific forward primer (10 mol/L), 0.4 uL of 
miScript universal primer (10 pmol/L), and 1 uL of PCR 
template (CDNA). The amplification protocol was as follows: 
initial denaturation and enzyme activation for 3 min at 95°C, 
followed by 40 cycles for 5 sec at 95°C, and for 20 sec at 60°C. 
Fluorescence was read at the end of each round of 
amplification. Each assay was performed in duplicate. The 
relative expression levels of one piRNA in the ovary and testis 
were determined using the (1+E7)°/(1+ER)“™ method, and 
different piRNAs in the same tissue were determined using the 
2°) method (Pfaffl, 2001). The 5S small ribosomal RNA of P. 
olivaceus was used as the internal control, and 
piRNA_0002790 was used as the calibrator. The primers of 
each piRNA are listed in Table 1, with the most conserved 


sequence of each PRC selected as the primer. The groups 
were compared via one-way analysis of variance (ANOVA) and 
Tukey's test to identify statistically distinct groups. Significant 
differences were accepted using P<0.05. 


RESULTS 


Histological observation and identification of piRNAs of gonads 
In the ovary of P. olivaceus, only mature oocytes or those close 
to maturity were revealed. The  previtellogenic oocytes 
displayed strongly basophilic ooplasms, while mature oocytes 
were strongly eosinophilic. In the testis, tubules were filled with 
spermatozoids, while the numbers of spermatogonia and 
spermatocytes were reduced (Figure 1). 
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Figure 1 Histological features of the ovary and testis of one-year-old P. olivaceus 
A: Ovary. 1, Stage Ill; 2, Stage VI; 3, Stage V. B: Testis. b, blood vessel; f: flagella; arrow indicates interstitial cells. 


Combined RNA samples from different ovary and testis 
samples were used for sequencing. A total of 24 399 395 and 
26 417 330 raw reads were obtained from the ovary and testis 
libraries, respectively. After filtering low-quality sequences, 
empty adaptors, and single-read sequences, 12 246 743 
(50.19%) and 13 260 806 (50.19%) clean reads were used for 
further sequence analyses. After assembly and BLAST 
searching in the piRNA database, 5 181 and 3 915 unique 
piRNA sequences showing high homology with zebrafish were 
obtained from the two libraries. All unique piRNAs were 
grouped into 289 piRNA clusters (PRCs), where all unique 
piRNAs in the same PRC shared a conserved sequence with 
varying lengths. The remaining nucleotides were similar to other 
types of RNAs, including ncRNA, tRNA, miRNA, rRNA, snRNA, 
and snoRNA. The size distribution of small RNAs for 
sequencing was similar between the two libraries. The majority 
of these small RNAs increased from 18 nt to 30 nt, and lengths 
ranged from 23 nt to 28 nt (Figure 2A). 


Expression patterns of PRCs in P. olivaceus gonads 

The relative abundance of various PRCs was determined using 
the small RNA deep sequencing approach by calculating 
sequencing frequency. In the ovary and testis libraries, the 
identified PRCs exhibited a broad range of expression levels. 
The most abundant PRC in the mature gonads was 
piR_0000766 (which displayed more than 12 000 reads in the 
ovary and testis), followed by piR_0034879 (10 766 reads). In 


contrast, 87 piRNAs displayed less than 10 reads. These 
results indicate that expression varied significantly among 
different piRNA families. 


Differential expression of PRCs in the ovary and testis 

On the basis of changes in the relative piRNA abundance 
between the two libraries, we found that log? € 5559^vY) of the 
expressed PRCs varied from -7.65 to 4.89. Approximately 
45.32% and 15.57% of isolated PRCs were preferentially 
expressed in the testis and ovary, respectively. Results showed 
a total of 289 PRCs, including 38 ovary-specific, 45 ovary-bias, 
24 testis-specific and 131 testis-bias PRCs (Figure 2B). Seven 
PRCs (piR_0048450, piR_0006070, piR_0019030, 
piR_0015958, piR_0055443, piR_0056629, and piR_0050867) 
with different reads obtained from sequencing results were 
selected and quantified using real-time RT-PCR to validate the 
results obtained from the Illumina deep sequencing data. These 
PRCs were also homologues, according to the piBASE 
database (http: //www.regulatoryrna.org/database/piRNA/). 
Three individual experiments using different RNA pools were 
performed to generate reliable and repeatable results. Results 
showed that PRCs, including piRNR 0048450, piRNA 0006070, 
piRNA 0019030, and  piRNA 0015958, exhibited high 
expression levels in both the testis and ovary, though were 
more highly expressed in the testis (P«0.05). Three other PRCs 
(piRNA 0055443, piRNA_ 0056629, and piRNA_ 0050867) 
were primarily detected in the ovary (P«0.05). These results are 
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consistent with the sequencing data (Figure 3). 
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Figure 2 Length range and differential expression of piRNA in 
the gonads of P. olivaceus 

A: Length distribution of sequenced small RNAs from the ovary and 
testis of P. olivaceus. X-axis denotes the piRNA length of 26-31 
nucleotides; Y-axis represents the percentage of piRNA. B: Scatter plot 
expression level of piRNA in testis versus ovary. Genes are listed in 
Supplemental Table 2. 


DISCUSSION 


As the largest group of vertebrates, teleosts exhibit remarkably 
varied sexuality (Godwin, 2010). Genetic, environmental, and 
social factors can affect sexual development in a large number 
of fish species at different developmental stages (Heule et al., 
2014). Paralichthys olivaceus, which is a common flatfish raised 
in China, Korea, and Japan, exhibits a temperature-sensitive 
mechanism of gender determination and differentiation. Artificial 
breeding and selection is performed to facilitate the production 
of P. olivaceus, and sex control is regarded as an important 
technique for fish production (Zhu et al., 2006). Although the 
sexual development of P. olivaceus needs to be properly 
understood, information concerning the involved mechanism 
remains insufficient and unclear. One-year-old Japanese 
flounder are suitable for supply to the market. The gonad is 
mature and easy to collect, observe and study in regards to 
gene expression. 

Three processes take place during sexual development, that 
is, sex determination, differentiation, and maintenance. The first 
process determines whether the bipotential primordium will 
develop into a testis or ovary. The second process involves the 
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Figure 3 Relationship between relative expression levels of 
piRNA validated by qRT-PCR and fold changes derived using 

small RNA sequencing 

Expression levels of piRNA measured with qRT-PCR are depicted with 
column charts, and triangles or circles represent the fold changes in 
sequencing. Error bars represent one standard deviation of three 
different biological replicates. Statistical significance is reported for each 
piRNA. F: female; M: male; *: P«0.05. 


actual development of the testes or ovaries from the 
undifferentiated gonads. The final process ensures the 
reproductive function of the ovary and testis (Freedberg & 
Taylor, 2007; Hayes, 1998; Spencer & Janzen, 2014). The 
protein-coding genes involved in the regulatory network of 
sexual development described so far originate from mammalian 
or Teleostei studies and can be divided into different functional 
groups. The dmy gene has been identified in medaka and plays 
a key role in gender development (Matsuda et al., 2002); other 
genes, including dmrt3 (Li et al., 2008), dmrt6 (Zhang et al., 
2014), gsdf (Shibata et al., 2010), 118-hydroxylase cyp11b 
(Schiffer et al., 2015), and 5a-reductase 1, 2, and 3 (Aquila et 
al., 2015), are also involved in sexual development. In addition, 
some genes, such as dmyldmrt1 in medaka, possess different 
biological functions in different strains (Matsuda et al., 2002). 
Plasticity is the main characteristic of sexual development, 
especially for fish. Studies on noncoding genes, including 
piRNA, microRNA, and IncRNA, can provide insights into 
sexual identification in fish (Mishima et al., 2008; Gu et al., 2014; 
Jing et al, 2014; Xiao et al, 2014) Gender-specific 
organization of piRNA has been discovered in Drosophila (Li et 
al., 2009; Malone et al., 2009), as have differences between 
male and female piRNAs in zebrafish (Houwing et al., 2007). 
Through microfluidic chip analysis, piRNAs have been found to 
exhibit quantitative differences, ranging from 1.58 to 64 019.26. 
A subset comprising more than 2096 of piRNA exhibited evident 
differential expression between the testis and ovary (Houwing 
et al., 2007). We also observed gender-based differential and 
specific organization of piRNA in P. olivaceus. Results suggest 


that the gender-related organization of piRNA might be 
conserved in vertebrates and flies, and might possess multiple 
functions in organisms. Piwi-antisense piRNA complexes 
mediate cleavage of sense piRNA precursors and transposon 
(or protein-coding) transcripts, which silence transposon and 
gene expression at the post-transcriptional level (Weick & 
Miska, 2014). Therefore, study on piRNA is important for 
illuminating the function of sex development-related protein- 
coding genes. 

piRNAs are involved in the sexual development of many 
vertebrates, particularly mammals and teleosts (Aravin et al., 
2007; Girard et al., 2006; Houwing et al., 2007). However, the 
related mechanism remains only partially understood because 
of the lack of information on piRNA expression patterns in the 
gonads of non-model organisms. Furthermore, identification of 
piRNAs is difficult because they lack conservative secondary 
structure motifs, such as the ones found in miRNA (Zhang et al., 
2010). Thus, we examined all potential piRNAs in both male 
and female gonads through deep sequencing technology, 
followed by BLAST searching the piBANK database with the 
zebrafish piRNA library as a reference. With a high frequency of 
24-28 nt RNA, the deep sequencing results suggested that 
piRNAs are present in the male and female gonads of P 
olivaceus. These findings are consistent with piRNA length 
distribution in zebrafish. Mammals, including humans, mice, 
and rats, have a major distribution interval of 27—31 nt (http: 
//pirnabank.ibab.ac.in/stats.html). The crystal structure of the 
MID domain is obtained from a Piwi/Argonaute protein, which 
plays a central role in RNA silencing process as essential 
catalytic components of the RNA-induced silencing complex 
(RISC). The RISC complex is responsible for the gene silencing 
phenomenon known as RNA interference. MID docking 
experiments showed its capability to specifically recognize the 
5'-uridine (1U-bias) of piRNAs (Cora et al., 2014). In this study, 
the ratios of A, T (U), C, and G in the first nucleotides of the 5 
865 unique piRNAs of P olivaceus were 23.2496, 30.7596, 
22.9596, and 23.0596, respectively, which also characterized 
1U-bias of the molecule. The expression levels of PRC in P. 
olivaceus presented a broad range. For example, more than 
several thousand reads of piR 0000766 were found in both 
male and female piRNA libraries, whereas specific infrequent 
PRCs showed several reads only We suggest that the 
imbalance in the expression levels of different piRNAs might 
correspond to the different abundance levels of their target 
genes. Hence, the relationship between the expression levels 
of PRCs and their target genes should be considered in future 
studies. 
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ABSTRACT 


Sphaeroma terebrans, a wood-boring isopoda, is 
distributed worldwide in tropical and subtropical 
mangroves. The taxonomy of S. terebrans is usually 
based on morphological characteristics, with its 
molecular identification still poorly understood. The 
number of teeth on the uropodal exopod and the 
length of the propodus of the seventh pereopod are 
considered as the major morphological characteristics 
in S. terebrans, which can cause difficulty in regards 
to accurate identification. In this study, we identified 
S. terebrans via molecular and morphological data. 
Furthermore, the validity of the mitochondrial 
cytochrome c oxidase subunit | (COI) gene as a 
DNA barcode for the identification of genus 
Sphaeroma, including species S. fterebrans, S. 
retrolaeve, and S. serratum, was examined. The 
mitochondrial COI gene sequences of all specimens 
were sequenced and analysed. The interspecific 
Kimura 2-parameter distances were higher than 
intraspecific distances | and no intraspecific- 
interspecific distance overlaps were observed. In 
addition, genetic distance and nucleotide diversity (77) 
exhibited no differences within S. terebrans. Our 
results revealed that the mitochondrial COI gene can 
serve as a valid DNA barcode for the identification of 
S. terebrans. Furthermore, the number of teeth on 
the uropodal exopod and the length of the propodus 
of the seventh pereopod were found to be unreliable 
taxonomic characteristics for S. terebrans. 


Keywords: Sphaeroma terebrans; DNA barcode; 
COI gene; Molecular identification 


INTRODUCTION 


Mangroves are biologically and globally important ecosystems 
(Giri et al., 2011). Their aerial roots provide an important 
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substrate in which many species of animals live and reproduce 
(Nagelkerken et al., 2008). Sphaeroma terebrans, a wood- 
boring isopoda, is found worldwide in tropical and subtropical 
mangroves (Estevez, 1978), where it preferentially burrows into 
the aerial roots for shelter and reproductive habitat (Harrison & 
Holdich, 1984; John, 1970). In recent years, substantial S. 
terebrans outbreaks have seriously affected mangrove stands 
in China, especially in Hainan island (Fan et al., 2014). 

The effects of S. terebrans on mangroves have been studied 
by many researchers (Estevez & Simon, 1975; Estevez, 1978; 
Jones & Icely 1981; Kensley & Schotte, 1999; Perry, 1988; 
Rehm & Humm, 1973); however, the taxonomic standards of S. 
terebrans remain poorly understood. Due to some minor 
morphological differences, including the number and 
arrangement of the tubercles on the pereonite, the structure of 
the pereopod, and the presence of tubercles furnished with 
bristle-like hairs on the abdomen, S. terebrans was previously 
named as S. vastator (Bate, 1866) and S. destructor 
(Richardson, 1897). Based on morphological identification, 
Estevez & Simon (1975) concluded that S. vastator and S. 
destructor were synonyms of S. terebrans. 

The classic use of morphological characteristics for species 
delimitation can result in under- or over-estimation of 
biodiversity due to factors such as phenotypic plasticity 
(Knowlton, 1993). DNA barcode, which can supplement 
taxonomic datasets in the process of species delimitation 
(Schindel & Miller, 2005), is a practical tool that can be used for 
the identification of various species within a known taxonomic 
framework and for linking different biological life stages of the 
same species (Feng et al, 2011; Puillandre et al., 2009; 
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Schindel & Miller, 2005). The mitochondrial COI gene has been 
proposed as a universal barcode, and has been successfully 
applied in the identification of Portunidae, fish, bivalve molluscs, 
and hoverflies (Blair et al., 2006; Hebert et al., 2003; Ma et al., 
2012; Persis et al., 2009; Stahls et al., 2009). The COI gene 
sequences of S. terebrans have been analysed in America and 
Africa (Baratti et al., 2005, 2011), with results suggesting that 
cosmopolitan S. terebrans is comprised of more than one 
species. Therefore, its taxonomic status needs to be revaluated. 

The aim of the present study was to provide a reliable and 
valid way to delimit S. terebrans. In this study, the validity of the 
mitochondrial COI gene as a DNA barcode marker for the 
identification of three species of Sphaeroma, namely, S. 
terebrans, S. retrolaeve, and S. serratum, was examined. To 
detect if there was any cryptic species in S. terebrans, the COI 
gene sequences of individuals with morphological differences 
were analysed. Our study should be useful in the identification 
of the genus Sphaeroma and for further research on S. 
trerebrans. 


MATERIALS AND METHODS 


Sampling and scoring of morphological characteristics 

The S. terebrans specimens were collected from three localities 
in China (Figure 1). Prior to DNA extraction, all specimens were 
examined under an anatomical lens and assigned to groups 
according to comparison with previous morphological 
descriptions (Harrison & Holdich, 1984). The morphological 
differences of S. terebrans were then photographed by a 
TM3030Plus tabletop microscope. The S. terebrans individuals 
were sorted according to the following morphological 
characteristics: the number and arrangement of tubercles on 
the pereonite, number of teeth on the uropodal exopod, shape 
of the pleotelson, setae distribution, and length of the second 
and seventh pereopods. These are considered to be major 
characteristics for the diagnosis of S. terebrans within 





Figure 1 Sample collection sites of S. terebrans 
HK: Haikou, Hainan, WC: Wenchang, Hainan, BH: Beihai, Guangxi, 
ZJ: Zhanjiang, Guangdong 
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Sphaeroma (Harrison & Holdich, 1984). The S. retrolaeve 
specimens were collected from Hainan and Beihai mangroves. 
All samples were preserved in 95% alcohol. 


DNA extraction, PCR amplification, and sequencing 

The genomic DNA of S. terebrans and S. retrolaeve were 
obtained from the pereopods. DNA extractions were performed 
using a TaKaRa MiniBEST Universal Genomic DNA Extraction 
Kit Ver.5.0 following the manufacturer's protocols. The primers 
mtd10 5-TTGATTTTTTGGTCATCCAGAAGT-3' (Roehrdan. 
1993) and Florence 5-CCTAAAAAATGTTGAGGGAA-3' were 
used for amplification of the mitochondrial COI gene (Baratti et 
al., 2005). We followed PCR protocols as per Baratti et al. 
(2005). The PCR products were electrophoresed using 196 
agarose gel and sequenced by Shanghai Majorbio Bio-Pharm 
Technology Co., Ltd. 


Data analysis 

All sequences were aligned using ClustalW (Thompson et al., 
1997). Interspecific and intraspecific sequence divergences 
were calculated using the Kimura 2-parameter (K2P) model 
with the pairwise deletion option in MEGA 5.0 (Kimura, 1980). 
Haplotypes were identified and analysed using DNA SP version 
4.1 (Librado & Rozas, 2009). Nucleotide diversity (rr) and 
haplotype diversity(h) were calculated according to Nei (1987) 
using DNA SP version 4.1 (Rozas & Rozas, 1999). Based on 
the K2P model, neighbor joining (NJ) and maximum likelihood 
(ML) trees were constructed using MEGA 5.0 (Kimura, 1980; 
Tamura et al., 2011), with the Cymodoce fuscina voucher from 
the NCBI (GenBank Accession No. KJ410468) used as an 
outgroup. Node supports for the two approaches (NJ and ML) 
were inferred with bootstrap analysis (1 000 replicates). 


RESULTS 


Prior to DNA extraction, we assigned S. terebrans specimens 
into A1-A5 and B morphotypes (Figure 2). The number of teeth 
on the uropodal exopod and the length of the propodus of the 
seventh pereopod varied within S. terebrans, which were 
assigned into seven groups (Table 1). 

Partially aligned COI sequences 498 bp in length were 
obtained from 70 S. terebrans individuals and 10 S. retrolaeve 
individuals. The COI sequences of S. serratum and C. fuscina 
voucher were downloaded from the NCBI. Details of these 
sequences are shown in Table 2. There were fifteen haplotypes 
for S. terebrans and two for S. retrolaeve. Haplotype sequences 
were deposited in the NCBI under accession numbers 
KU558703-KU558719. 

All haplotype sequences were aligned and edited, and no 
insertion or deletion sites were found in any of the sequences. 
The intraspecific distances in S. terebrans ranged from 0.001 to 
0.013 (Table 3). The maximum interspecific distance (1.394) 
was between S. serratum and C. fuscina voucher, while the 
minimum interspecific distance (0.24) was between S. serratum 
and S. retrolaeve. No overlaps between interspecific and 
intraspecific distances were found, suggesting the existence of 
a distinct barcoding gap. The NJ phylogenetic tree is shown in 


Figure 3. Distinct clusters corresponding to species were found 
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Figure 2 Diagnostic morphological characteristics of S. terebrans 





with high bootstrap support. 
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A1-A5 are uropodal exopods with different numbers of teeth. B is the seventh pereopod with different propodus length. 


The S. terebrans individuals were sorted into seven groups 
according to their morphological traits, and partial COI 
sequences of S. ferebrans were aligned and compiled. The 
intraspecific distances ranged from 0.001 to 0.003 within the SS, 
SW, WW, WL, and LL groups. The intraspecific distance 
between PL and PS was 0.001 (Table 4). The mean haplotype 
diversity (h) was 0.55596, and ranged from 0.20096 (PL group) 
to 0.866% (WW group) (Table 5). The highest nucleotide 
diversity (77) was found in the WW group (0.004), while the 
lowest was found in the PL group (0.000) (Table 5). Results 
suggested that there were no mitochondrial genetic variations 
within S. terebrans. 


Phylogenetic analysis of genus Sphaeroma was performed 
using NJ and ML methods, which yielded similar results. The NJ 
tree revealed that the three species of Sphaeroma and one 
species of Cymodoce formed monophyletic clusters (Figure 3). 
The nearest relationship was observed between S. retrolaeve 
and S. serratum, while the most distant relationship was found 
between S. ferebrans and C. fuscina voucher. 


DISCUSSION 
The rapid and effective identification of closely related wood- 


borer Sphaeroma species is important for the research and 
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Table 1 List of sampling localities and morphological differences of S. terebrans 


Teeth on uropodal exopod 








Codes Number Location Propodus 
Left Right 
SS1-5 5 Beihai, Guangxi 3 3 Short 
SS6-10 5 Wenchang, Hainan 3 3 Short 
SW1-5 5 Zhanjiang, Guangdong 3 4 Short 
SW6-10 5 Haikou, Hainan 3 4 Short 
WWh1-10 10 Haikou, Hainan 4 4 Short 
WL1-10 10 Wenchang, Hainan 4 5 Short 
LL1-5 5 Beihai, Guangxi 5 5 Short 
LL6-10 5 Haikou, Hainan 5 5 Short 
PS1-10 10 Wenchang, Hainan 4 4 Short 
PL1-10 10 Wenchang, Hainan 4 4 Long 


The numbers of teeth on the left and right exopod of S. terebrans are 3 and 3, respectively, for SS; 3 and 4 for SW; 4 and 4 for WW; 4 and 5 for WL; 
and 5 and 5 for LL. PL means that the propodus of the seventh pereopod of S. ferebrans is long, whereas PS is short. The same in the following table. 


Table 2 List of COI sequences, GenBank accession numbers, and geographic sources of samples 





Species n GenBank accession no. Source 

S. terebrans 15 KU558705-KU558719 Hainan, China 
S. retrolaeve 2 KU558703, KU558704 Guangxi, China 
S. serratum 1 GU130256 Germany 

C. fuscina voucher 1 KJ410468 Germany 


n indicates the number of sequences analyzed. 


Table 3 Pairwise genetic distances (Kimura 2-parameter) of three Sphaeroma species and one Cymodoce species based on COI sequences 


SR-H1 SR-H2 ST-H1 ST-H2 ST-H3 ST-H4 ST-H5 ST-H6 ST-H7 ST-H8 ST-H9 ST-H10 ST-H11 ST-H12 ST-H13 ST-H14 ST-H15 SS CF 





SR-H1 

SR-H2 0.004 

ST-H1 0.254 0.248 

ST-H2 0.251 0.245 0.006 

ST-H3 0.260 0.254 0.008 0.006 

ST-H4 0.251 0.245 0.008 0.006 0.008 

ST-H5 0.257 0.251 0.006 0.004 0.002 0.006 

ST-H6 0.254 0.248 0.004 0.002 0.004 0.004 0.002 

ST-H7 0.257 0.251 0.011 0.008 0.011 0.011 0.008 0.006 

ST-H8 0.257 0.251 0.006 0.004 0.006 0.006 0.004 0.002 0.008 

ST-H9 0.251 0.245 0.011 0.004 0.006 0.011 0.004 0.006 0.013 0.008 

ST-H100.251 0.251 0.006 0.004 0.006 0.006 0.004 0.002 0.008 0.004 0.008 

ST-H110.254 0.248 0.006 0.004 0.006 0.006 0.004 0.002 0.008 0.004 0.008 0.004 

ST-H120.254 0.248 0.006 0.004 0.006 0.006 0.004 0.002 0.004 0.004 0.008 0.004 0.004 

ST-H13 0.254 0.248 0.004 0.002 0.004 0.004 0.002 0.001 0.006 0.002 0.006 0.002 0.002 0.002 

ST-H14 0.254 0.248 0.004 0.002 0.004 0.004 0.002 0.001 0.006 0.002 0.006 0.002 0.002 0.002 0.001 

ST-H15 0.257 0.251 0.008 0.006 0.008 0.008 0.006 0.004 0.011 0.006 0.011 0.006 0.006 0.006 0.004 0.004 

SS 0.241 0.240 0.243 0.246 0.249 0.243 0.246 0.243 0.252 0.246 0.246 0.246 0.240 0.246 0.243 0.243 0.243 

CF 1.290 1.270 1.099 1.131 1.111 1.135 1.123 1.123 1.135 1.123 1.144 1.135 1.114 1.135 1.123 1.123 1.131 1.394 
H1-H15: Haplotype 1-15, SR: S. retrolaeve, ST: S. terebrans, SS: S. serratum, CF: C. fuscina voucher. 
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KU558711 
KU558716 
KU558714 
KU558719 
KU558710 
KU558708 
KU558717 
KU558718 
KU558706 
KU558713 
KU558707 
KU558709 
KU558712 
KU558715 
KU558705 


S. terebrans 




















GUI130256 | S. serratum 
KU558703 
100! KU558704 
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S. retrolaeve 





re 


0.1 


Figure 3 
haplotypes of three species of Sphaeroma and one species of 
Cymodoce 


Neighbor-joining phylogenetic tree of individual 


Table 4 Intraspecific genetic distances (Kimura 2-parameter) of S. 


terebrans with morphological differences based on COI sequences 


SS SW ww WL LL PL PS 





SS 

SW 0.002 

WW 0.003 0.003 

WL 0.002 0.002 0.003 

LL 0.001 0.002 0.003 0.002 

PL 0.001 0.001 0.002 0.002 0.001 

PS 0.001 0.002 0.003 0.002 0.001 0.001 


Table 5 Number of haplotypes, haplotype diversity(h), and nucleotide 
diversity(7) of different groups 





Group No. Haplotype Nucleonde 
haplotypes diversity (h) diversity (77) 
SS 3 0.377 0.001 
Sw 5 0.666 0.002 
WW 7 0.866 0.004 
WL 6 0.777 0.002 
LL 3 0.600 0.001 
PL 2 0.200 0.000 
PS 3 0.416 0.001 


restoration of eroded mangroves. Identification of S. terebrans 
based on morphological characteristics alone is weak and, to 
some extent, ambiguous. Based on morphological characteristics, 
Some individuals of S. terebrans were previously named S. 
vastator (Bate, 1866) and S. destructor (Richardson, 1897). In 


this study, clear evidence was provided for the identification of S. 


terebrans individuals, which exhibited differences in morphological 


characteristics. The validity of using the mitochondrial COI gene 
sequence as a DNA barcode for the identification of genus 
Sphaeroma was examined, and included three Sphaeroma 
species, namely, S. terebrans, S. retrolaeve and S. serratum, 
with C. fuscina voucher (Sphaeromatidae) used as an outgroup. 
A distinct barcoding gap was found between the intraspecific 
and interspecific distances in each species. The NJ 
phylogenetic tree consisted of four distinct clusters, each 
containing individuals from one species only. These results 
indicate that the partial mitochondrial COI gene is an effective 
DNA barcode for the identification of the genus Sphaeroma. 

Individuals of S. terebrans had different numbers of teeth on 
the uropodal exopod and different lengths of the propodus of 
the seventh pereopod. These individuals were sorted into 
seven groups, with each group containing 10 individuals. The 
genetic distance and nucleotide divergence showed no 
variation among the different groups. Therefore, these results 
revealed that the COI gene sequences of individuals with 
morphological differences were almost no difference. Although 
Harrison & Holdich (1984) determined that the propodus of the 
seventh pereopod of subadult males is relatively short, Our 
investigations showed that the length of the pereopodal 
propodus in S. terebrans was not necessarily linked with gender. 
Previous research concluded that cosmopolitan S. terebrans 
was comprised of more than one species (Baratti et al., 2011, 
2005), but morphological taxonomic details of S. terebrans were 
not mentioned. In our research, specimens in China were 
carefully checked according to morphological characteristics 
and were assigned into different groups, with molecular 
methods used for further identification. This combination of 
morphological taxonomy and molecular divergence should 
provide results of greater reliable. 


CONCLUSIONS 


In this study, the mitochondrial COI gene was found to be an 
effective DNA barcode for the identification of Sphaeroma 
species, whereas the number of teeth on the uropodal exopod 
and the length of the propodus of the seventh pereopod were 
found to be invalid taxonomic characteristics. The phylogenetic 
relationships determined in this study will be of use for studying 
the species composition of Sphaeroma in eroded mangroves in 
China and for establishing a good foundation for the restoration 
of mangrove ecosystems. 
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Z Visiting Professor, Hainan Medical University, Haikou Hainan 571101, China 


DEAR EDITOR: 


We read the publication by Liu (2016) on authorship misconduct 
with a great interest. In fact, the authorship misconduct is no 
uncommon and has to be managed. Liu (2016) noted that 
“Proper authorship embodies honesty, integrity, fairness and 
transparency, which surely are the very essence of any 
scientific pursuit.” Due to the present requirement on academic 
publication to support academic rank, several forms of 
authorship misconduct can be seen. Of interest, some 
publications might have more than twenty authors which is an 
extremely high number. To give the authorship as a gift to the 
senior staff or professor is not uncommon practice in many 
developing countries. Sometimes, a worse situation, senior staff 
might force the junior to put his/her name in the author’s list. How 
to manage these situations is the big question in medical ethics. 
Some senior one(s) might tell that the work has already been 
thought by him/her/them and the junior has to give the authorship 
to him/her/them. To avoid dispute, the junior usually allow and 
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permit as the senior(s) desire(s) (Rajasekaran et al., 2015). 

In addition, the more serious situation can also be seen. 
Sometimes, the submitted article contains plagiarisms or frauds, 
many authors in the author list try to deny and present 
his/her/their innocence. Many senior staffs might point to the 
junior author and claim that the junior author use his/her/their 
names without permission. So the system to support the good 
authorship clarification practice is needed (Rajasekaran et al., 
2015; Mandal et al., 2015). 
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Appointment of Dr. Wai-Yee Chan as the Associate 
Editor-in-Chief of Zoological Research 


We are very pleased to announce that we have invited Dr. Wai-Yee Chan, Professor and Director of 
School of Biomedical Sciences, The Chinese University of Hong Kong, to serve as the Associate Editor- 
in-Chief for Zoological Research (ZR), effective 1 September, 2016. 

Dr. Chan’s research interests include developmental genomics and epigenomics of germ cell 
development, germinal stem cell biology, functional genomics of developing male gonads, molecular 
genetics of human reproductive and endocrine disorders, non-coding RNA, antisense transcript, 
microRNA in normal and abnormal development, and genetics and genomics of foetal development. 
Many of his findings have been published in high-impact peer-reviewed journals such as Science, Proc 
Natl Acad Sci USA, Nucl Acid Res, Hum Mol Genet. 

i In his new position, with his impressive academic achievements and well-recognized expertise, Dr. 
À Chan will certainly help promote the international impact of ZR. 
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